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Layered  plates  fall  into  several  categories  where  the  preliminary  assumptions 
as  well  as  the  methods  of  analysis  are  different,  thus  leading  to  multitudinous  of 
specialized  theories.  Starting  from  the  three-dimensional  theory  of  non-linear 
elasticity,  and  using  a  displacement  approach,  much  like  in  the  Reissner-Mindlin  type 
theories,  a  general  layered  plate  theory  is  developed.  It  has  the  advantage  of 
producing  a  complete  set  of  fundamental  equations  consistent  with  various  stages  of 
linearization  in  the  general  strain-displacement  relations.  The  notion  of  stresses  in 
the  reference  state  is  employed  in  the  fundamental  equations,  and  a  variational 
procedure  is  extended  to  derive  the  equations  of  motion  and  the  boundary  conditions 
of  layered  plates  subjected  to  large  displacements  and  large  angles  of  rotation. 

The  present  theory  incorporates  the  effects  of  transverse  shear  and  transverse 
normal  strains  as  well  as  rotatory  inertia,  with  different  thickness,  material  densities 
and  material  constants  in  each  layer  of  the  plate.  All  stress  components  in  each  layer 
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are  considered,  and  the  continuity  conditions  between  the  contiguous  layers  are 
maintained  for  all  tractions.  Each  layer  is  assumed  to  be  anisotropic,  having  elastic 
symmetry  with  respect  to  its  middle  plane.  The  assumed  displacement  field  is 
piecewise  linear  in  the  thickness  coordinate  in  all  displacement  components  and 
fulfills  the  geometric  continuity  conditions  between  the  contiguous  layers. 

Theoretical  development  of  the  mathematical  form  of  the  nonlinear  equations 
of  motion  and  the  boundary  conditions  is  presented  for  the  case  of  three-layered 
plate.  In  demonstrating  the  generality  of  the  theory,  several  theories  available  in  the 
literatures  are  derived  from  the  present  as  special  cases.  A  numerical  investigation 
is  conducted  to  correlate  the  results  of  the  linear  theory  with  a  simplified  version  of 
the  present  theory;  and  the  effects  of  compressible  layers  on  a  multi-layered  plate 
instability  are  examined. 
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CHAPTER  I 
INTRODUCTION 


1.1  Introduction  to  The  Theory  of  Plates 

The  classical  plate  theory  came  as  a  result  of  the  significant  treatment  of 
plates  in  the  1800s  principally  by  Navier,  Kirchhoff,  Levy,  Lagrange  and  Cauchy. 
Based  on  the  fundamental  assumptions  known  as  the  Kirchhoff  hypotheses,  the 
theory  is  limited  to  the  small  deflection  and  bending  of  isotropic,  homogeneous, 
elastic,  thin  plates  [1].  In  the  case  of  thin  plates  with  large  deflection,  Kirchhoff  s 
second  assumption  (as  defined  by  Love  [1])  is  satisfied  only  if  the  plate  is  bent  into 
a  developable  surface  [2].  Otherwise,  bending  of  a  plate  is  accompanied  by  strain 
in  the  middle  plane.  The  corresponding  stresses  must  be  taken  into  consideration 
in  deriving  the  differential  equation  of  plates.  In  this  way  we  obtain  non-linear 
equations,  and  the  solution  of  the  problem  becomes  much  more  complicated  [2]. 
Two  coupled  non-hnear,  partial  differential  equations  governing  the  large  deflection 
problem  were  first  introduced  by  von  Karman  in  1910  [3].  The  approximate  theories 
of  thin  plates  become  unreliable  in  the  case  of  plates  of  considerable  thickness.  In 
such  a  case,  the  thick-plate  theory  which  considers  the  problem  of  plates  as  a  three- 
dimensional  problem  of  elasticity,  should  be  applied.  The  stress  analysis  becomes, 
consequently,  more  involved.  Thus  Kirchhoffs  first  assumption  (also  known  as 
Kirchhoff  s  hypothesis)  becomes  no  longer  valid.  However,  it  was  first  recognized 
that  transverse  shear  strains  in  the  plate  can  not  be  neglected  when  plate  vibration 


2 

was  studied  [4].  In  1951,  Mindlin  [5]  developed  his  plate  theory  to  include  transverse 
shear  in  a  similar  manner  to  Timoshenko's  beam  theory  of  1921  [6]. 

A  trend  toward  a  theory  of  plates  which  do  not  possess  an  isotropic  behavior 
later  gained  larger  interest.  Multi-layered  plates  is  one  example  of  such  behavior. 
The  simplest  approach,  In  the  beginning,  to  these  types  of  plate  structures  was  to 
assume  a  single  layer  orthotropic  plate.  In  general,  flexural  and  extensional  motions 
are  coupled  in  a  layered  plate.  Composites  were  also  analyzed  using  similar  theories 
[7].  However,  sandwich  plates,  which  were  very  popular  thirty  some  years  ago, 
attracted  more  attention;  hence,  several  sophisticated  approaches  led  to  a  level  of 
sandwich  plate  theories  in  a  class  by  itself. 

Laminated  composites  have  been  for  a  long  time  and  are  still  being  analyzed 
as  anisotropic  plates  of  the  classical  type,  anisotropy  being  the  major  additional 
feature  that  was  accounted  for  [7].  For  symmetrical  laminated  plates  flexure  and 
extension  are  uncoupled.  If  the  bending-twisting  stiffness  coupling  terms  in  the  plate 
equation  further  vanish,  the  symmetric  laminates  are  referred  to  as  specially 
orthotropic  and  were  usually  treated  as  orthotropic  plates,  also  of  the  classical  type. 
But  since,  recently,  the  transverse  shear  effect  has  been  shown  to  have  a  great 
influence  on  the  overall  behavior,  whether  dynamics  and  vibration  or  bending  and 
buckling  analysis  are  concerned,  the  analysis  of  laminated  composite  plates  has  taken 
new  directions  from  the  classical  theory  of  laminates  [7].  In  contrast,  the  shear  effect 
was  shown  to  be  important  in  the  analysis  of  sandwich  plates  a  long  time  ago; 
sandwich  plate  theories  have  developed  very  well  and  several  higher  order  theories 
have  been  introduced.   Only  recently,  upon  this   development,   the    connection   of 
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laminated  composites  to  layered  plate  analysis  in  general  has  been  examined  and 

exploited   [7].     Thus,  the  much  advanced  theories  in  sandwich  plates  became 

attractive  tools  to  the  analysis  of  layered  plates  in  general  and  laminated  composites 

in  particular.    In  the  following  section  we  will  examine  this  trend  from  its  early 

stages. 

1.2  A  Prospective  of  the  Theories  of  Multi-Layered  Plates 

The  expression  "sandwich  plate"  designates  a  composite  plate  consisting  of  two 
very  thin  layers  of  high-strength  "face"  material,  between  which  a  thick  layer  of  ultra- 
lightweight  "core"  is  sandwiched.  The  early  sandwich  plate  structures  had  identical 
thickness  for  the  faces  which  preserved  symmetry  about  the  midplane.  By  the  late 
1940s  Hoff  and  Mautner  [8]  had  derived  differential  equations  and  boundary 
conditions  for  the  bending  and  buckling  of  sandwich  plates.  They  have  reflected  in 
their  derivation  the  general  understanding  of  the  behavior  of  sandwich  plates  at  that 
time.  The  moment  of  inertia  of  the  cross  section  of  the  sandwich  plate  is  large 
because  of  the  comparatively  great  distance  between  the  two  faces,  and  thus  the 
buckling  formulas  derived  from  classical  plate  theory  give  buckling  stresses  that,  as 
a  rule,  far  exceed  the  yield  stress  of  the  face  material.  Since  the  modulus  of  elasticity 
of  the  core  in  the  plane  of  the  plate  is  of  the  order  of  magnitude  of  one-thousandth 
of  that  of  the  faces,  the  normal  stresses  in  the  core  are  of  little  importance  in 
resisting  bending  moments  even  though  the  usual  ratio  of  face  thickness  to  core 
thickness  is  between  one-tenth  and  one-hundredth.  The  lightweight  core  permits 
unusually  large  shearing  deformations,  and  performs  a  task  of  transmitting  shear 
forces.    Hence  shearing  deformations  must  not  be  disregarded  in  the  analysis  of 
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sandwich  plates.    Moreover,  relative  displacements  of  the  two  faces  are  possible 

because  of  the  small  extensional  rigidity  of  some  of  the  core  materials  used. 

The  theory  of  Hoff  and  Mautner,  with  the  aid  of  the  principle  of  virtual 

displacements,  is  developed  from  a  consideration  of  these  strain  energy  portions 

stored  during  deformations,  the  strain  energy  caused  in  the  faces  by  extension  and 

bending  and  that  caused  in  the  core  by  shear  and  by  extension  perpendicular  to  the 

plane  of  the  faces.  Later  in  1950,  Hoff  [9]  modified  the  theory  and  added  the  strain 

energy  of  in-plane  shear  in  the  faces,  but  continued  to  disregard  that  of  extension 

parallel  to  the  faces  in  the  core.    He  justified  the  assumption  that  the  in-plane 

stresses  in  the  core  contribute  only  negligible  amounts  to  the  total  strain  energy,  by 

taking  the  elastic  and  shear  moduli  of  the  core  to  be  small  as  compared  with  those 

of  the  faces.    Also,  normal  strains  in  the  core  are  disregarded.    Furthermore,  the 

strain  energy  stored  in  the  faces  because  of  shear  perpendicular  to  the  faces  is 

neglected  -  this  is  permissible  when  the  ratio  of  the  length  or  width  of  the  plate  to 

the  thickness  of  a  face  is  always  large.  The  displacement  field  was  described  by  three 

functions  corresponding  to  the  coordinate  system  directions.  The  in-plane  deflections 

are  of  the  faces  and  they  are  equal  and  opposite  for  each  face,  while  the  vertical 

deflection  for  the  entire  plate  takes  place  through  shearing  of  the  core.  This  vertical 

deflection  does  not  cause  force  resultants  (corresponding  to  membrane  stresses)  in 

the  individual  faces  but  it  gives  rise  to  bending  and  twisting  moments  in  them 

because  of  the  non-vanishing  bending  and  torsional  rigidity  of  the  faces.    Small 

deflection  elasticity  is  assumed,  and  the  material  of  each  layer  is  taken  to  be 

isotropic.  One  last  note  to  Hoft^s  work  is  that  he  acknowledged  that  the  contribution 

of  the  core  to  resisting  bending  cannot,  as  a  rule,  be  neglected. 
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In  1948,  E.  Reissner  introduced  a  large  deflection  theory  which  extends  Foppl- 

von  Karman  theory  of  large  deflection  of  ordinary  plates  to  sandwich  plates  [10].  His 

equations  permit  the  analysis  of  the  effect  of  transverse  shear  stress  deformation  and 

transverse  normal  stress  deformation  in  the  core  on  the  overall  behavior  of  the  plate. 

Reissner's  results  led  him  to  an  important  conclusion  at  the  time.   It  states  that  the 

range  of  deflections  for  which  the  linear  "small  deflection"  theory  is  adequate 

decreases  in  accordance  with  a  simple  explicit  formula  as  the  core  is  made  softer 

relative  to  the  faces.   This  was  also  encountered  by  Hoff  s  experimental  results. 

Reissner  extended  the  two  assumptions  that  the  thickness  of  each  face  tf  is 

small  compared  with  that  of  the  core  t,  and  that  the  value  of  the  elastic  constants  E^, 

Gf  for  the  face  layers  are  large  compared  with  the  values  of  the  elastic  constants  E^, 

G(,  for  the  core  layer,  to  further  assume  that  the  products  tfEf,  tfOf  are  large 

compared  with  the  values  of  t^E^.,  and  tfi^.  On  the  basis  of  the  first  assumption,  he 

assumes  that  the  stresses  in  the  faces  parallel  to  their  planes  are  distributed 

uniformly  over  the  thickness  of  the  face  layer.  On  the  basis  of  the  third  assumption, 

he  neglects  the  face-parallel  stresses  in  the  core  layer  and  their  effect  on  the 

deformation  of  the  composite  plate.  As  a  result,  the  theory  treats  the  sandwich  plate 

as  a  combination  of  two  plates  without  bending  stiffness  (the  face  layers),  and  of  a 

third  plate  (the  core  layer)  offering  resistance  only  to  transverse  shear  stresses  and 

transverse  normal  stresses.  The  displacement  field  is  described  by  six  displacement 

functions,  three  components  in  each  face  layer  corresponding  vdth  the  coordinate 

system.    An  interesting  note  to  be  made  is  that  Reissner  had  opted  to  follow  an 

approach  parallel  to  the  method  of  elasticity  in  deriving  his  system  of  differential 

equations.     That  is  he  started  by  integrating  the  equations  of  equilibrium  and 
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consequently  he  had  to  define  several  other  terms  in  his  theory.   Then,  through  an 

elaborate  algebraic  manipulation  the  stress-strain  and  the  strain-displacement 
relations  are  included,  and  the  result  is  a  set  of  several  algebraic  relations  along  with 
a  system  of  non-linear  partial  differential  equations  for  three  displacement  functions 
and  an  additional  three  functions  defining  the  changes  of  slopes. 

In  the  meantime,  various  others  had  been  working  on  the  problem  of 
sandwich  plates.  However,  various  basic  assumptions  were  made,  such  as  (1)  the 
core  undergoes  shear  deformation  only;  (2)  the  bending  rigidity  of  the  core  is 
neglected;  (3)  some  of  the  shear  and  normal  stresses  which  occur  in  the  core  are  also 
neglected;  and  (4)  the  faces  are  assumed  to  be  membranes. 

In  1952,  Eringen  presented  a  theory  of  bending  and  buckling  of  sandwich 
plates  which  became  one  of  the  important  works  in  the  area  and  one  which  would 
last  for  years  [11].  To  the  present  day  it  is  still  being  used  as  an  important  reference. 
The  interest  in  this  work  came  mainly  from  the  fact  that  the  above  mentioned 
restrictions  are  removed.  Consequently,  the  theory  is  more  general  in  that  it 
considers  all  six  components  of  the  stress  tensor  in  the  core,  face  layers  having 
bending  rigidity,  and  it  encompasses  two  types  of  bending  and  buckling,  namely: 
overall  bending  and  buckling  and  bending  and  buckling  with  flattening 
(compressibility)  of  the  core. 

The  simplifying  assumptions  are  within  the  scope  of  the  usual  small 
deformation  theory,  i.e.  (1)  faces  are  thin  as  compared  to  the  core,  but  are  not 
membranes-thus,  for  the  bending  of  the  faces  the  Bernoulli-Navier  hypothesis  is  valid 
(plane  sections  perpendicular  to  the  median  plane  of  the  plate  remain  plane  during 
deformation);  (2)  the  displacements  are  linear  functions  of  the  distance  from  the 
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median  plane  of  the  plate;  and  (3)  the  linear  theory  of  elasticity  is  assumed.   Since 

the  faces  are  taken  to  be  thin  and  symmetrical,  he  uses  one  set  of  displacement 

components  for  the  entire  plate  while  assuming  linearity  in  distribution  up  to  the 

median  plane  of  the  faces.  In  other  words,  assumption  (2)  is  the  result  of  expanding 

the  displacement  components  into  a  power  series  of  the  transverse  coordinates  and 

taking  into  account  only  the  terms  which  are  linear  in  the  transverse  coordinate. 

In  the  analysis,  the  total  potential  of  the  sandwich  plate  is  expressed  in  terms 
of  the  deflection  components,  and  thus  with  the  aid  of  the  principle  of  the  minimum 
of  the  total  potential  energy,  four  partial  differential  equations  are  obtained.  The 
usual  methods  of  the  variational  calculus  are  used  in  order  to  obtain  the  extremals 
of  the  total  potential.  Eringen  next  solved  an  example  problem  using  Navier  method 
and  Fourier  series  type  solution. 

It  is  clear  that  the  simplicity  of  Eringen's  theory,  while  not  compromising  some 
of  the  generality  of  the  problem,  is  due  mainly  to  its  linearity.  One  aspect  of  the 
linear  theory  is  assuming  linear  distribution  for  the  deflection  components  through 
the  entire  plate,  whence,  limiting  the  application  of  the  theory  to  sandwich  plates 
with  very  thin  face  layers. 

Y.  Y.  Yu  generalized  Eringen's  concept  by  assuming  thick  faces  and  using 
different  displacement  functions  for  each  layer,  though  this  was  done  indirectly.  Yu 
had  started  from  three-dimensional  elasticity  and  followed  the  procedures  of  small 
deformation  theory  as  defined  by  Novozhilov  [12].  Yu's  contribution  to  the  problem 
of  sandwich  plates  in  particular  and  that  of  layered  plates  in  general  spans  over  thirty 
years  and  includes  a  long  list  of  published  literatures.  It  started  in  1959  when  he 
introduced  a  new  theory  of  sandwich  plates  [13].     By  assuming  the  transverse 
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displacement  to  be  uniform  across  the  plate  thickness,  and  the  displacements  in  the 

plane  of  the  plate  to  be  of  linear  variations,  with  the  slope  in  the  faces  taken  to  be 
different  from  that  in  the  core,  it  was  possible  to  include  in  the  theory  the  effects  of 
transverse  shear  deformations  in  the  core  and  faces,  the  rotatory  and  translator/ 
inertias  of  the  core  and  faces,  the  flexural  rigidity  of  the  core,  and  the  flexural  and 
extensional  rigidities  of  the  faces.  Hence,  although  the  sandwich  as  a  whole  is  under 
flexure,  each  of  the  two  face  layers  is  under  combined  flexure  and  extension.  Since 
no  restrictions  were  imposed  on  the  magnitudes  of  the  ratios  of  the  thicknesses  and 
elastic  constants  between  the  core  and  face  layers,  the  results  were  therefore  in 
general  applicable  to  any  symmetrically  arranged  three-layered  plate.  Continuity  of 
tractions  as  well  as  displacements  is  maintained  at  the  interfaces  between  adjacent 
layers.  In  Yu's  early  work,  the  non-linear  equations  of  motion  of  sandwich  plates 
together  with  the  appropriate  boundary  conditions  were  derived  from  the  variational 
equation  of  motion  of  the  non-linear  theory  of  elasticity.  In  his  later  work  [14,15], 
Hamilton's  principle  in  dynamics,  which  is  the  counterpart  to  the  minimum  total 
potential  energy  principle,  was  applied  in  the  derivations. 

In  all  of  the  foregoing  investigations  the  core  and  facings  of  the  sandwich 
plates  were  assumed  to  be  isotropic.  Alwan  [16]  in  1964  introduced  an  analysis  that 
adopted  the  same  assumptions  as  Reissner  [10]  but  the  core  was  taken  as  an 
orthotropic  honeycomb-type  structure. 

While  Yu  may  be  considered  one  of  the  first  to  adopt  a  continuous  piecewise 
linear  displacement  distribution,  the  concept  of  piecewise  displacement  distribution 
has  become  an  obvious  representation  in  the  theory  of  layered  plates.  Another 
breakthrough  came  in  the  mid-  to  late  1960s  when  Ebcioglu  and  Habip  introduced 
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the  non-linear  equations  of  motion  of  plates  and  shells  in  the  reference  state  [17,18]. 

The  theory  of  the  reference  state  involves  the  notion  of  stress  measured  per  unit  area 

of  the  undeformed  body,  taken  as  a  reference,  in  contradistinction  from  the 

conventional  representation  of  stress  measured  per  unit  area  of  the  deformed  body. 

The  distinction  is  of  special  significance  for  the  correct  interpretation  of  the 

"geometrically  non-hnear,"  "finite-"  or  "large  deflection"  theories  of  plates  and  shells, 

where  this  point  is  often  overlooked.    The  non-linear  equations  of  motion  in  their 

theory  had  been  obtained  by  integrating  the  corresponding  three-dimensional  stress 

equations  of  motion  through  the  thickness  of  the  undeformed  body.  The  concept  of 

a  reference   state  makes  it  possible  to  properly  identify  the  origin  of  various 

additional  terms  that  appear  in  the  field  equations  of  non-linear  theories  of  plates 

and  shells  in  comparison  with  those  of  linear  theories.   Alternatively,  Habip  gave  a 

new  derivation  of  these  equations  of  a  non-linear  theory  of  elastic,  anisotropic  and 

heterogeneous  plates  and  shells  by  means  of  the  modified  Hellinger-Reissner 

variational  theorem  of  three-dimensional  continuum  dynamics  [19,20]. 

However,  the  further  extension  of  these  considerations  to  the  case  of  elastic 

sandwich  plates  is  due  to  Ebcioglu  [21,22].  He  considers  a  sandwich  plate  taking  into 

account  all  the  stress  components  in  each  layer.  As  in  other  works  by  Yu  [13,14],  this 

theory  incorporates  the  effects  of  transverse  shear  and  transverse  normal  strains  as 

well  as  rotatory  inertia,  with  different  material  constants  in  each  layer  of  the 

sandwich  panel.    It  is  assumed  that  the  facing  and  the  core  are  anisotropic,  each 

having  elastic  symmetry  with  respect  to  its  middle  plane.  An  added  consideration  to 

the  theory  was  including  the  effect  of  steady  thermal  gradients  in  the  stress-strain 

relations.     In  a  system  of  convected  general  curvilinear  coordinates,  the  stress 
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equations  of  motion  are  incorporated  into  the  variational  integrals  which  stem  from 

the  application  of  Hamilton's  principle.  In  the  displacement  field,  it  is  no  longer 
assumed  that  the  transverse  displacement  component  is  uniform  across  the  plate 
thickness,  and  the  angle  of  rotations  are  assumed  to  be  different  in  each  layer.  The 
continuity  of  tractions  and  displacements  at  interfaces  are  preserved. 

Ebcioglu  extended  his  geometrically  non-linear  theory  of  sandwich  plates  to 
the  range  of  material  non-linearity  [23]  (physical  non-linearity,  in  the  sense  of 
Novozhilov  [12]).  Simplified  non-linear  strain-displacement  relations  are  used.  For 
the  elastic  case,  each  layer  is  of  a  different  thickness  and  of  a  different  anisotropic 
material  having  one  plane  of  elastic  symmetry.  Transverse  shear,  rotatory  inertia  and 
thermal  effects  are  included  in  each  layer.  However,  the  influence  of  transverse 
normal  strain  on  the  deflection  is  neglected.  For  the  plastic  case,  the  stress-strain 
relations  of  the  Henky-von  Mises  deformation  theory  of  plasticity  are  used  including 
the  effect  of  temperature  and  compressibility.  For  the  latter  case,  only  isotropic 
materials  are  assumed. 

Vasek  and  Schmidt  [24]  introduced  a  theory  for  the  non-linear  bending  and 
buckling  of  multi-sandwich  plates.  Their  plate  consisted  of  N  stiff  layers  and  N-1 
weak  layers.  The  assumptions  made  were  as  follows:  (1)  the  linearly  elastic, 
homogeneous,  isotropic  stiff  layers  possess  bending  stiffness  and  deform  in 
accordance  with  Kirchhoffs  hypothesis  for  thin  plates;  (2)  the  linearly  elastic, 
homogeneous,  orthotropic  weak  layers  can  transmit  the  transverse  normal  and 
shearing  stresses  but  not  the  in-plane  stresses;  (3)  the  layers  behave  non-linearly  in 
a  sense  of  Foppl-von  Karman;  (4)  the  thicknesses  and  the  materials  may  be  different 
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from  layer  to  layer;  and  (5)  the  layer-to-layer  bonds  are  strong  enough  so  that  under 
all  loadings  no  bond  failure  will  occur. 

In  a  series  of  studies,  including  one  on  laminated  and  sandwich  plates  in  1970 
[25],  Pagano  constructed  three-dimensional  elasticity  solutions  and  compared  them 
to  the  analogous  results  in  classical  laminated  plate  theory.  He  concluded  that  the 
laminated  plate  theory  (LPT)  leads  to  a  very  poor  description  of  laminate  response 
at  low  span-to-depth  ratios,  but  converges  to  the  exact  solution  as  this  ratio  increases. 
This  convergence,  he  observed,  is  more  rapid  for  the  stress  components  than  plate 
deflection. 

Based  on  the  same  conclusion,  laminated  plate  theory  based  on  the  Kirchhoff 
hypothesis  is  inaccurate  for  determining  gross  plate  response  and  internal  stresses  of 
thick  composites  and  sandwich  type  laminates.  Whitney  in  1972  [26]  introduced  a 
procedure  which  is  an  extension  of  the  LPT  to  include  the  effect  of  transverse  shear 
deformation.  In  1979,  Bert  [27]  stated  that  thickness-shear  deformation  is  important 
even  in  single-layer  panels  of  composite  material  due  to  their  very  low  ratio  of  shear 
modulus  to  Young's  modulus  as  compared  to  isotropic  materials. 

A  departure  from  the  approach  of  Pagano  and  Whitney  of  considering  all 
layers  as  one  equivalent  single  anisotropic  layer  was  presented  by  Di  Sciuva  in  1987 
[28].  He  assumed  the  displacement  field  to  be  piecewise  linear  in  the  in-plane 
components  and  uniform  in  the  transverse  component. 

Before  we  depart  from  our  review,  we  must  mention  one  last  important 
contribution  to  the  laminated  plates  analysis.  It  is  the  1984  Reddy's  higher  order 
theory  for  laminated  composite  plates  [29].  As  was  stated  earlier,  laminated 
composite  plates  were  until  recently  analyzed  as  anisotropic  plates.  First-order  shear 
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deformation  plate   theories  emerged,  whicli   consider  linear  distribution  of  the 

displacement  of  the  laminate  and  the  linear  strain-displacement  relations  of  the  small 
deformation  elasticity.  The  stress-strain  relations  are  then  incorporated  in  the  layers, 
and  hence  the  laminate  is  defferentiated  from  an  anisotropic  plate. 

Reddy  attempted  to  improve  the  accuracy  of  the  prediction  of  stresses  and 
displacements  in  the  laminate  while  not  increasing  the  number  of  dependant 
variables  and  satisfying  the  boundary  conditions  at  the  surfaces.  To  maintain  the 
satisfaction  of  the  condition  that  the  transverse  shear  stresses  vanish  on  the  surfaces 
and  are  non-zero  elsewhere,  a  parabohc  distribution  of  the  transverse  shear  strain  is 
required.  Reddy  then  required  the  use  of  a  displacement  field  in  which  the  in-plane 
displacements  are  expanded  as  cubic  functions  of  the  thickness  coordinate  and  the 
transverse  deflection  is  constant  through  the  plate  thickness.  He  justified  the  uniform 
distribution  of  the  transverse  deflection  by  comparing  the  in-plane  and  transverse 
normal  stresses.  The  number  of  generalized  coordinates  is  reduced  by  setting  to  zero 
the  transverse  shearing  stresses  at  the  top  and  bottom  surfaces.  The  linear  strain- 
displacement  relations  are  also  used,  while  the  plate  is  assumed  incompressible,  i.e., 
transverse  normal  strain  is  ignored.  The  stress-strain  relations  in  each  layer  possess 
a  plane  of  elastic  symmetry.  The  stress  continuity  across  each  layer  interface  was  not 
imposed.  The  theory  is  then  characterized,  according  to  Reddy,  as  a  simple  two- 
dimensional  theory  of  plates  that  accurately  describes  the  global  behavior  of  the 
laminated  plates  and  seems  to  be  a  compromise  between  accuracy  and  ease  of 
analysis. 

Several  scholars  attempted  to  improve  on  Reddy's  theory;  most  notable  are 
Pandaya  and  Kant  [30],  and  Librescu,  Kudeir  and  Frederick  [31].    Several  refined 
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anisotropic  composite  laminated  theories  substantiated  on  the  basis  of  different  initial 

assumptions  have  been  considered  and  compared  by  Librescu  and  Reddy;  and  it  was 

shown  analytically   [32]   and  numerically  [33]   that  they  represent  but  different 

formulations  of  a  single  theory,  designated  as  the  moderately  thick  plate  theory. 

The  higher-order  shear  deformation  theories,  as  well  as  the  first-order  theory 
will  not  fulfill  the  continuity  conditions  for  the  transverse  shearing  stresses  at  the 
interfaces  [28]. 

On  the  numerical  analysis  application  to  the  theories  of  layered  plates,  the 
most  notable  work  is  the  finite  element  method  implementation  into  the  solution  of 
the  problem  of  layered  plates.  The  earliest  attempt  was  introducing  the  triangular 
element  for  multi-layer  sandwich  plates  by  Khatua  and  Cheung  [34].  Recently, 
Rajagopal,  Singh  and  Rao  [35]  introduced  a  quadratic  isoparametric  element  having 
five  degrees  of  freedom  at  each  node  for  the  non-linear  analysis  of  sandwich  plates. 
Limited  to  linear  conditions,  a  composite  finite  element  analysis  requiring  only 
continuity  of  the  field  variables,  thus  permitting  the  use  of  simple  C°  elements  as 
opposed  to  the  requirements  for  the  more  complicated  finite  element  analysis,  was 
presented  by  El-Hawary  and  Herrmann  [36]. 

The  preceding  review  of  numerical  analyses  of  layered  plates  by  no  means  is 
representative  of  the  computational  work  in  this  area.  Clearly,  our  main  concern  is 
the  theoretical  and  analytical  advancements  in  the  field;  and  the  review  of  these 
advancements  is  unequivocally  complete  as  of  the  date  of  this  study.  To  summarize, 
layered  plates  fall  into  three  categories,  sandwiches,  multi-sandwiches  and  laminated, 
as  illustrated  in  Figure  1.1.  The  preliminary  assumptions  as  well  as  the  methods  of 
analysis  in  each  of  these  categories  are  different,  thus  leading  to   multitudinous   of 
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LAYERED  PLATES 


SANDWICHES 


MULTI-SANDWICHES 


STIFF-SOFT-STIFF 


LAHINATED 


ALTERNATING 
STIFF-SOFT 
LAYERS 


STIFF  LAYERS  -  VERY  THIN 
WEAK  LAYERS  -  THICK 


OVERALL  BEHAVIOR 
THIN  LAYERS 


I 


IDENTICAL  THICKNESS 


/  PRESENT 
I  THEORY 


NO  THICKNESS  OR  MATERIAL  LIMITATIONS 


Figure  1.1:  Categories  of  Layered  Plates. 

specialized  theories.  The  present  theory  is  a  general  theory  in  the  sense  that  it  is  not 
limited  to  any  one  of  these  categories,  rather  it  is  applicable  to  layered  plates  in 
general. 

1.3  Principles  of  the  Present  Theory 

Starting  from  the  three-dimensional  theory  of  non-linear  elasticity,  and  using 
a  displacement  approach,  much  like  in  the  Reissner-Mindlin  type  theories,  we  will 
develop  our  three-dimensional  layered  plate  theory.  The  present  work  can  be 
considered  as  an  extension  of  the  theory  of  sandwich  plates  formulated  by  Ebcioglu 
[21,22],  and  as  a  refinement  of  the  contribution  to  the  theory  of  sandwich  plates  by 
Eringen  [11]  and  by  Yu  [13,14].   The  notion  of  stresses  in  the  reference  state,  used 
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by  Habip  [19]  and  Ebcioglu  [21,22],  is  employed  in  the  fundamental  equations.  Also, 

a  variational  procedure  used  by  Yu  [15]  is  extended  to  derive  the  equations  of 

motion    and    the    boundary    conditions    of    layered    plates    subjected    to    large 

displacements  and  large  angles  of  rotation. 

The  present  theory  incorporates  the  effects  of  transverse  shear  and  transverse 

normal  strains  as  well  as  rotatory  inertia,  with  different  material  densities  and 

material  constants  in  each  layer  of  the  plate.  All  stress  components  in  each  layer  will 

be  considered,  and  no  assumption  underlying  the  significance  of  a  stress  component 

is  implemented  anywhere  in  the  derivation.  The  continuity  conditions  between  the 

contiguous  layers  is  maintained  for  the  stress  vectors  (tractions).   In  addition,  each 

layer  may  be  of  a  different  thickness,  and  no  a  priori  limitations  are  imposed  upon 

the  displacement  functions  which  define  collectively  the  displacement  field,  which, 

otherwise,  would  produce  appreciable  constraints  during  the  deformation  of  the  plate 

as  discussed  by  Ebcioglu  [21],  especially  for  the  case  of  symmetrical  flattening 

discussed  by  Eringen  [11].    Each  layer  is  assumed  to  be  anisotropic,  having  elastic 

symmetry  with  respect  to  its  middle  plane.    The  assumed  displacement  field  is 

piecewise  linear  in  the  thickness  coordinate  in  the  in-plane  as  well  as  the  transverse 

components  and  fulfills  the  geometric  continuity  conditions  between  the  contiguous 

layers;  furthermore,  it  takes  into  account  the  distortion  of  the  deformed  normal. 


CHAPTER  II 
GENERAL  THEORY 


2.1  Mathematical  Preliminaries 

Let  every  point  of  the  continuous  three-dimensional  body,  called  briefly  the 
body  B^,  be  at  rest,  at  time  t=^t^  relative  to  a  fixed  rectangular  Cartesian  system  of 
axes  X,  (see  Appendk  A).  The  position  vector  of  a  typical  point  P^  of  the  body  5„ 
referred  to  the  origin  is  given  by 

r=^A-'A-  (2.1) 

where  i^.  are  unit  vectors  along  the  fixed  axes. 

We  suppose  the  body  B^  is  deformed  so  that  at  time  t  a  typical  point  F„  has 

moved  to  P.   The  position  vector  of  P  referred  to  the  same  origin  is 

R  =  yk  h  (2.2) 

The  position  vector  of  F  relative  to  P^  is  denoted  by  V  and  is  called  the  displacement 
vector.   Thus 

V  =  R-r  =  {y,-x,)i,  (2.3) 

We  assume  that  each  point  P,  at  time  t,  is  related  to  its  original  position  F„  at  time 
t-t^hy  the  equations 

y.  =  yi  (xi  ,  ^2  ,xs,t)  (2.4a) 

Xi  =  x^  (yi ,  yj .  ^i  >  0  (2.4b) 

where  >',  andx,  are  single-valued  and  continuously  differentiable  with  respect  to  each 
of  their  variables  as  many  times  as  may  be  required.   If  this  deformation  is  to  be 
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possible  in  a  real  material  then 

[|ij>0  (2.5) 

The  general  theory  of  the  present  work  is  based  on  the  derivation  in  terms  of 
the  reference  state  of  the  plate.  Thus,  the  rectangular  Cartesian  system  of  axes 
defined  in  the  undeformed  body  at  time  t=tg  are  used  as  the  base  system  for  all 
vectors  and  tensorial  quantities. 

The  base  vector  g,  and  metric  tensor  g,y  may  be  defined  for  the  coordinate 
system  x,  in  the  body  jS„,  so  that 

g,  =  r,;  =  i/  (2.6a) 

g,j  =  g,  •  g;    =  ^ij  (2.6b) 

Similarly,  base  vector  G,-  and  metric  tensor  G,j  may  be  defined  for  the  coordinate 
system  jc,  in  the  body  B  at  time  t.   Thus 

G,  =  R,,  (2.7a) 

««  =  «.•«-  =  5- a^  P.7b) 

where  the  metric  tensor  G^  is  called  in  some  literature  the  Green's  Deformation 
tensor  when  it  is  referenced  to  the  coordinate  system  in  the  undeformed  body,  as  is 
the  case  here,  instead  of  a  general  curvilinear  system  d-  as  an  example. 
The  symmetric  strain  tensor,  Ytp  is  defined  by  the  equation 

y,  =  V2  (  G^.  -  g,^ )  (2.8) 

From  (2.3)  we  see  that        G,-  =  R,^  =  r,.  +  V„- 
Hence,  using  (2.6),  (2.7),  and  (2.8), 

/,■  =  V2  (  g,  .  V,^.  +  gj .  V,  +  V,,  .  V,^.  )  (2.9) 

The  displacement  vector  V  may  be  expressed  in  terms  of  the  base  vectors  oiB^.  Thus 


18 

V  -  V,„  g,„  V,,  =  V„„  g,,,  (2.10) 

By  introducing  (2.10)  into  (2.9)  and  using  (2.6),  the  strain  tensor  y,j  is  expressed  in 
terms  of  the  displacement  components  and  becomes  as 

y.   =  1  (V.  .+  V.  .  +  V  .  y    )  (2.11) 

where  V,.  are  the  components  of  displacement  referred  to  the  axes  .v,  in  the 
undeformed  body. 

2.2  Strain-Displacement  Relations 

The  general  non-linear  strain-displacement  relations  for  large  elastic 
deformations  are  given  in  tensor  form  in  (2.11)  which  are  identical  to  those  given  in 
most  literature  on  elasticity  and  non-linear  continuum  mechanics,  e.g.[37].  The  terms 
of  (2.11)  denote  three-dimensional  space  functions.  However,  in  the  study  of  plates 
and  shells,  it  is  more  convenient  to  refer  to  surface  functions  (without  losing  any  of 
the  transverse  and  out  of  plane  terms.)   Thus  we  write 

ya3  =  \iK^,-yy.-yra  ^r^^^v^  ^3^3)  (2.12) 

Y33-^(2K3,3.K,3K,3.F3,3K3,3) 

Expressions  (2.12)  are  identical  to  (2.11). 

In  the  theory  of  large  deflections  of  thin  plates,  a  simplified  non-linear  strain- 
displacement  relations  have  been  used  extensively,  i.e.  [3], 

ya,    -i_iK^p-yp,a-^V.    ^vp)  (2-13) 

These  are  based  primarily  on  the  following  assumptions:  (1)  plane  sections  before 
deformation  which  are  perpendicular  to  the  middle  surface  of  the  plate,  remain  plane 
and  perpendicular  to  the  middle  surface  after  deformation;  and  (2)  the  distance  of 
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every  point  of  the  plate  from  the  middle  surface   remains  unchanged  by  the 

deformation  [3].    However,  in  developing  a  general  theory  of  large  deflections  of 

multi-layered  plates,  the  transverse  shear  and  normal  strains,  (y^.-,  y^-,  and  y..)  can 

not  be  ignored.   Consequently,  the  above  assumptions  become  invalid. 

A  simplified  non-linear  form  of  (2.12)  has  been  used  in  the  theories  of 

sandwich  plates  [21], 

V      =   ^  (V     +V     +V      V    ) 

Y33    ~    ^3'3 

where  the  assumption  that  the  maximum  deflection  is  assumed  to  be  of  the  order  of 
magnitude  of  the  thickness  is  removed.  The  relations  (2.14)  will  be  used  in  the 
present  theory.  The  in-plane  relations  are  non-linear,  while  the  transverse  shear  and 
normal  strain-displacement  relations  are  hnear,  and  this  is  by  no  means  unjustifiable 
since  these  components  are  defined  in  each  layer  rather  than  throughout  the  entire 
plate  thickness. 

A  more  general  form  than  (2.14)  has  been  used  in  [18]  when  considering  a 
large  deflection  theory  of  anisotropic  plates, 

2.3  Stress-Strain  Relations 

The  stress  vector  ^t,  per  unit  area  of  the  undeformed  body,  associated  with  a 
surface  in  the  deformed  body,  whose  unit  normal  in  its  undeformed  position  is  ^n,  is 

«t  -  S,j  „«,■  Gj  (2.16a) 

o"  =  o^i^i  (2- 16b) 
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The  stress  tensor  5^  (second  Piola-Kirchhoff  stress  tensor)  is  measured  per  unit  area 

of  the  undeformed  body  while  defining  the  state  of  stress  in  the  deformed  body. 

The  equations  of  motion  in  terms  of  5^  are 

i"  ''o   ■■  [S.r  (5,,.   +    Vj,r)h    +    Po  LPj)    =    9o  Uj\  (2.17) 

where,  F  =  /-,  g,,      and      f  =  ^  g, 

are,  respectively,  the  body  force  vector  and  acceleration  vector,  p^  is  the  density  of 
the  undeformed  body,  and  v^  is  the  volume  of  the  undeformed  body.  From  the 
moment  equations  of  motion,  it  follows 

5,.  =  5^,  (2.18) 

whence,  the  stress  tensor  5^  is  symmetric  [38]. 

For  an  elastic  body,  an  elastic  potential  or  strain  energy  function  W  exists, 
measured  per  unit  volume  of  the  undeformed  body,  which  depends  on  the  strain 
components  Vy  and  has  the  property 

5W^*  =  S,  5y,  (2.19) 

The  variations  <5Yy  are  subject  to  conditions 

^Yy  =  5y,,-  (2.20) 

so  that  for  a  compressible  body, 

5.  =  i(^.^)  (2.21) 

"     '    ^Y,       ay. 

If  the  stress  vector  ^t  be  referred  to  base  vectors  g,  in  the  undeformed  body, 

«t  =  t,j  on,  gj  (2.22) 

where,  /^,  =  5,,  (5^,  +  V^,,)  (2.23) 

is  another  stress  tensor  (first  Piola-Kirchhoff  stress  tensor)  measured  per  unit  area 

of  the  undeformed  body,  depicting  the  state  of  stress  in  the  deformed  body.    The 
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equations  of  motion  in  terms  of  /^  are 

tip.  +  Po  ioFj)  =  9o  iofj)  (2.24) 

where  tjj  is  not  symmetric  but  satisfies, 

t.,,,  g,„  •  Gj  =  tj,„  g„,  .  G,  (2.25) 

As  a  more  specific  set  of  stress-strain  relations  than  (2.21),  we  may  assume  the 
linear  version 

S.j  =  Cy„  Yr.  (2-26) 

The  stiffness  tensor  holds  the  following  symmetry  relations, 

For  a  medium  having  elastic  symmetry  with  respect  to  the  surface  Xj  =  constant 
(z  =  constant),  equations  (2.26)  reduce  to  the  following  form 


(2.28) 


Equations  (2.28)  express  a  similar  form  to  the  well-known  law  of  Hooke.  For 
infinitely  small  elongations  and  shears  (but  not  displacements  and  angles  of  rotations, 
which  are  in  no  way  restricted  by  the  preceding  assumption)  not  exceeding  the  limit 
of  proportionality,  the  stress-strain  relations  are  linear.  However,  the  above 
assumption  does  not  imply  a  linearized  version  of  the  strain-displacement  relations. 

2.4  Displacement  Field 

Referring  to  the  undeformed  plate,  we  will  define  a  rectangular  Cartesian 
coordinate  system  at  the  middle  surface  of  each  layer  as  shown  in  (Fig.  2.1).  The 
points  of  the  3-D  space  of  the  plate  will  be  referred  to  by  a  set  of  orthogonal 


s., 

^a/3YM    **'yM     "^    ^a/333    Y33 

Sa, 

-    ^a3Y3    Yy3 

^33 

^33yM    "^YM     "^    ^3333    Y33 
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coordinates  x„  where  x^,  (a  =  1,2)  and  x^ 
denotes  the  in-plane  coordinates  and 
the  normal  one  to  the  plane  x^-0, 
respectively.  The  flat  plate  consists  of 
A^  layers.  A  superscript  in  the  brackets 
"<>"  associated  with  a  field  quantity 
identifies  its  affihation  to  the  layer  of 
the  laminate  indicated  by  that  index. 
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igure  2.1:  Coordinates  System  in  Each 
Layer  with  Cross-Section  in  .Xj  -  X3  Plane. 


The  middle  layer  (core)  has  index  <0> 

and  layers  above  the  middle  plane  have  odd  indices  while  those  below  the  middle 
plane  have  even  indices.  Let  2//"^"^  be  the  thickness  of  layer  «  where  {n:n6[0,N-\]}. 
Thus  we  have 


^a         ^a 


<1> 


<2> 


^3    = 


=    ....    =  X^ 


<N-1> 


,;">  ^[h<o>+2  Y^    h^^-'^+h^"^],    oddn>l 


n/2 


^<n>  _  [/;<«> +25; /j<2^-^>+/i<">],      even  n>2 


r=2 


.  (2.29) 


We  now  assume  that  the  components  V^  of  the  displacement  vector  V  are 
given  in  each  layer  approximately  by 

where,  u^  and  ij;,  are  functions  ofx,y  and  t  only  and  have  the  dimensions  of  length 
and  slope  (dimensionless),  respectively.  The  displacement  components  (2.30)  contain 
six  unknown  displacement  functions  for  each  layer,  that  is  a  total  of  6*yV  unknowns 
in  the  plate.   These  can  be  reduced  in  number  to  3*(A^+ 1)  displacement  functions 
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by  implementing  the  conditions  of  displacement  continuity  at  interfaces, 

V-"\x^,x,=H<">+{-\rh<">-t) 

=  Vr'\x^,x,=H<"-'->  ^{-\r-'h<"-'->  ;t),    n=2X.,N-\ 
where, 


.    (2.31) 


.1-3  =  x;""-  +//<"> 


//<">  = 


Vz{n-\) 


+  [//<"> +2  Y.    /z<2'-i> +/,<"> ]^    odd  n 

/■=i 

-[/;■=">+ 2^ /;<-'-2>+/i<">],      even  n 


.        (2.32) 


The  conditions  (2.31)  are  applied  to  (2.30)  and  are  simplified,  hence 
t<;''^+[/;<°>t;°>+2  'Y.    /z^^-^^tr"^"'' Wi->i|;;"=-],  oddn 


U: 


;=1 
„/2 
..  ■^"■^      r;   <o>    I  <o>      T  V~v    I   <2;-''>    ■  <2r-2>      ,   <„>    ,  <«>, 

«,■       -["        Vi       +22^ /z"^-"  -^ij;.  +/z^"'^ij;.      ],  evenn 

r=2 


(2.33) 


Upon  substitution  of  (2.33)  into  (2.30),  we  obtain 


f/<0>  <0>  ,    <0> 

Vi         =   U-  +  X3II;, 


K""=< 


•/2(;i-l) 


"72 


(2.34) 


which  describe  a  linear  piecewise  displacement  representation  in  the  thickness 
coordinate  of  the  multi-layered  plate.  Since  the  strain-displacement  relations  are 
non-linear,  the  stress  quantities  are  quadratic  in  the  thickness  coordinate  (Reddy's 
dilemma  eliminated!) 

A  special  case  which  will  be  given  a  special  consideration  in  the  present  study 
is  a  three  layered  plate.  For  TV =3,  the  displacement  field  representation  is  given  by 
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v; 


u,  +  X3VI;,. 


T/<1>  <0> 


xr<2>  <0> 


<i>  ^i.  <i>i,i,<i> 


(2.35) 


where  superscripts  1,  0  and  2  indicate  upper  face  layer,  core  and  lower  face  layer, 
respectively,  as  shown  in  (Fig.  2.2).  The  displacement  components  in  (2.35)  are 
expressed  in  terms  of  the  12  displacement  functions 


/       <o>      ,  <o>      ,  <1>      ,  <2>  X 
(",  ,^,         .^,         >^/         ) 


(2.36) 


Higher     order     terms     in     the 

expansions    given    by    (2.35)    can    be 

included.  The  method  of  analysis  to  be 

presented  is  applicable  for  higher  order 

expansions      for      the      displacement 

components.  In  our  approximations,  by 

including  second  order  terms  for  the 

faces    in    addition    to    the    core    in 

transverse  direction,  we  have  added  the     Figure      2.2:      Rectangular      Cartesian 

Coordinates  in  Three-Layered  Plate. 

capability    to    allow   for   compressible 

faces  as  well  as  a  compressible  core.  In  that  respect,  the  present  displacement  field 
represent  an  extension  to  that  of  Ebcioglu  [22].  Additionally,  interface  continuity  is 
preserved  while  non-linear  distribution  through  the  thickness  of  the  panel  and  the 
differences  among  the  axial  in-plane  components  of  the  displacements  in  the  layers 
are  not  ignored. 
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2.5  Hamilton's  Principle--The  Variational  Integral 

In  dynamics,  the  counterpart  of  the  minimum  potential  energy  theorem  is 

Hamilton's  principle.    It  states;  The  time  integral  of  the  Lagrangian  function  over  a 

time  interval  tj  to  t2  is  stationary  for  the  "actual"  motion  with  respect  to  all  admissible 

virtual  displacements  which  vanish,  first,  at  instants  of  time  tj  and  r,  at  all  points  of  the 

body,  and,  second,  overs^,  where  the  displacements  are  prescribed,  throughout  the  entire 

time  interval  [39}.  Accordingly,  for  the  conservative  body  force  and  surface  traction, 

and  for  an  arbitrary  time  interval  t^  and  ?,,  the  Lagrangian  function  L  has  to  satisfy 


Ldt  =  0 


(2.37) 


where; 


L^K-  {U  +  A) 


U  = 


SA  = 


K  =  A 

7 


S,j  Yy  dv 


dV.  dV. 

'-dv 


J^SV.ds 


°  a    a 

The  volume  integrals  of  (2.38), 
(2.39),  and  (2.40)  represents, 
respectively,  the  strain  energy,  the 
variation  of  the  potential  energy  due  to 
the  body  forces  gb,  per  unit  mass  of  the 
undeformed  body,  and  the  kinetic 
energy,  and  must  be  extended  over  the 


(2.38) 
(2.39) 

(2.40) 


^igure      2.3:      Undeformed      body      in 
Reference  State. 
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volume  Vg  of  the  body  in  its  undeformed  state.    The  surface  integral  of  (2.39) 

represents  the  work  of  the  surface  forces  J^  (assumed  held  constant)  on  the  (virtual) 
variations  in  the  displacement  trajectory  [variations  from  the  actual  trajectory]. 

The  components  of  the  surface  force  J^  are  referred  to  base  vectors  in  the 
undeformed  body.  Let  „^  (underlined)  be  the  prescribed  value  of,/,  on  the  boundary. 
Then  the  stress  boundary  conditions  are 


oU  -  ot.  -  0 


(2.41) 


The  variational  integrals  of  the  stress  equations  of  motion  and  boundary 
conditions  consistent  with  strain-displacement  relations  (2.14)  are  obtained  from 
Hamilton's  principle.    We  can  write  (2.38)  in  the  following  form 


U  = 


i^apyap   +25„3Y„3+533Y33)^V 


(2.42) 


Based  on  (2.19),  property  of  the  strain  energy  function  iV,  it  can  be  shown  that 


6U  = 


(■^./3 -57,^+ 25,367,3+533 5 Y33)rfv^ 


(2.43) 


With  the  use  of  the  strain-displacement  relations  (2.14)  and  the  divergence  theorem, 

(2.44) 


'^.P^P  ^^    =^   ^^ap    Jh  ^^ 


we  put  (2.43)  into  the  following  form 


5t/  = 


\^o.p>p^y.^{^.pyy^p^y^^dv  -J(5,3,35k,  -5,3,, 5^3)^,. 

533,3  5K3rfv  -  \{^ap  ^K  ^  S,,  K3,,  5K3)  ^n^ds 


The  variation  (2.39)  is  written  as 


^^^^y^Jhds 


(2.45) 


SA   = 


PoioK^K  \b,8V,)dv  -\Ujv,  \i^sv.;)ds         (2.46) 
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And,  from  (2.40)  we  obtain 


Kdt  = 


(    p^V^6V^dv)dt 


rr  d 


[^{p^V^6V,-pXsK)]dv)dt 
ot 


[rp„F,5F,^v]-r( 

i\p^V.6V.dv)dt 


p^V^SV^dv)dt 


which  simply  takes  the  form, 


Kdt 


{PoK^^a- 9oy,&y,)dvdt 


(2.47) 


We  now  substitute   the   results  from  (2.45),   (2.46)  and   (2.47)   into  (2.37),  and 
Hamilton's  principle  in  a  form  consistent  with  plate  theory  notation  is 


'^^JU5,,,,-5,3,3.p,(A-^'J]5K„ 

-[i^a.yv.l,  -S^„^  -S33,3  +  p„(A  -  V,)]  6V,}dv 

I 

'-  4  -  (^^^^3>.  ^S,,ln^  -533  J,,]  6V,  ]ds  =  0 
where,  the  surface  forces  are 


Ja   =  ^aH   Jh   ^  '^a3  JH 


(2.48) 


(2.49) 


o^3  =  ('^a/)^3'.   ^^a-i^Jh   ^^^-iJh 

Here,  the  vector  {ji^  ,  /z^)  represents  the  unit  normal  associated  with  the 
undeformed  body  with  respect  to  the  base  vectors  in  the  undeformed  body.  An 
element  of  volume  dv  can  be  written  as 


dv  =  dz  dA 


(2.50) 
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where,  z  is  the  coordinate  in  the  transverse  direction  of  the  undeformed  body  where 

2  =  0  is  the  middle  plane  of  the  core,  and  dA  is  an  element  of  area. 

The  volume  integral  (2.48)  is  divided  into  the  number  of  the  layers  in  the 
laminate;  and  the  surface  integral  is  also  divided  to  correspond  to  the  number  of 
layers  in  addition  to  the  upper  and  lower  face.  Based  on  (2.50),  an  integration  in  the 
thickness  coordinate  may  be  performed.  This  procedure  leads  to  the  derivation  of 
the  equations  of  motion  and  boundary  conditions  as  will  be  demonstrated  in  the 
following  chapters. 

A  note  needs  to  be  made  here  regarding  the  expressions  in  (2.49).  These 
surface  force  expressions  are  extracted  from  the  derivation  of  the  variational  integral 
(2.48),  and  hence,  are  consistent  with  the  components  of  the  strain  field  given  in 
(2.14).  Since  the  relations  (2.14)  are  modified  non-linear  strain-displacement 
relations  and  do  not  constitute  a  complete  strain  tensor  as  may  be  defined  for  the 
general  non-linear  case,  traditional  tensorial  procedures  to  derive  the  surface  forces 
(i.e.,  Malvern  [40])  can  not  be  implemented  to  derive  (2.49). 


2.6  The  Constitutive  Equations 

Define  the  stress  resultants  in  each  layer  as  follows: 


iK7  M:i\Ki\Qr  jr.N^") 

.  <;■>     „  <r> 


<r>     <r>^    c<'>    c<i'>^<r>    c<r>^ 


(2.51) 


Define  body  force  resultants  with  respect  to  the  middle  planes  of  the  layers 


H,"' 


{Fr,Fr,Mr,Mr)=  r  pr^(x^^x^^x^^^^^^.^r^^^^=•)^^^'"  (-•-'-) 
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and  the  acceleration  resultants 


[.<r>    /-<r>         <r>         <r>N 


<r>  /^"T<r>  T'><r>  t><''> 


(/a       'A       ''"<.       ''"3       )  =  Po      (f^a       .^3       '^a       ^         '^3       ^         )^^ 


(2.53) 


The  edge  surface  force  resultants  with  respect  to  the  middle  plane  of  the  layer  are 


^-a      '^3       '-a      '-3       ''  "  ^ola      'o^       'o^a       ^         'o-3       ' 


')dz 


(2.54) 


where  j^^  ,  and  J_^  are  prescribed  traction  forces  on  the  surface  of  the  plate  as  were 
defined  in  the  proceeding  section. 

From  the  strain-displacement  relations  (2.14)  and  the  displacement  field 
representation  (2.34),  we  obtain  the  following  relations 


Ya/3       =oYa/3     +2  J^^     +2  .Ya^ 


<r>  <r> 


Ya3       =oYa3     +^  lY„3 


<r>  <r> 


'33 


=oY33 


where 


<o>     _ 

oYa^ 


<o> 
lYa^ 

<o> 
2Ya^ 


r  ,  <0>     ,  <0>       <0>  ,  <0>       <0>   ,  <0>T 

[ta,/3     +^^,a     +"3,a     ^3,^     +"3,/?     'I's.a     ] 

r      ,     <0>      ,     <0>   T 

[1'3,<.      ^3,^     ] 

^u  ^  i    r     ,    <0>  <0>  T 

0Ya3        =t[^c<        +"3,a      ] 

<o>  1    ,   ^^. 

lYa3        =7l|'3,a 

<o>  ,   <( 

0Y33         =  11^3 


(2.55) 


(2.56a) 
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<o>  ,  <o; 


<rn)  ■>  1,   <n>  <n>  <o>        <0>         i    ^n-->  /   ^    <0>  ,    <0>  <.o>,<.c'->      <.u^     .    ^u-'\ 

dY.T      =T[<r^"/r^"3;     "3,/?      ^^     "    i^a,P     ^^P,a     +"3,a      *3,^      +"3,/3      ^3,a     ) 
.   ^^->2      <o>       <o>       ,   <-;,ii/  .  <m>        .  <«!>  <o>    ,  <m>  <o>    ,  <in> 

+  /i       ^3,a     ^3,?     +/'  (ifa,/?      + 1*/?..*      +«3,a     ^3,^      +"3,^     "^Xa 

I    <o>    .  <o>    ,  <'">       ;   <o>    1  <o>    ,  <m>s    ,;,  </H>^.i<"' >,■<'">  1 
+/J  t3.a     ^3,^       +/i        ^3,/J      t3,a      ) +«  ^3,a       ^3,/!      J 


('"-l)/2 


E  [/' 


<2j-1>/  ,  <2'-l>        ,  <2r-l>        ,  <2r-l>      <o>  ^   .<>-!>  ,<o>  s 


/•=1 


</»>    ,   <2r-l> 


<;?;>    ,  <2j-1  > 


+  /'  (^3,a     ^3,^  +^3,/J     ^3.a  )+''  (^3,a      "^  XP  +  ^3,/3       ^3,<.  ) 

(m-l)/2  (/7,-l)/2 

+ 

(=1  i=l 

<»i>         1  r   .   <"'>        ,   <»!>  <o>    ,   <m>   ,       <o>    ,   <'">  ,  7,  <o>/,i  <o>  ,i,'="'>   ,   ,i,*=''>  ,|,'^"'^M 

lYc/J        =7[ta,/3      +^/3,c»      +^'3,«      I'S,/?      +"3,^      ^3,a      +''  i'^'^.a     1^3,/3      +  ^3./3     *3.<.      )J 

(m-l)/2 


y:  e  (/^^"-'^/'^^-^^t3:r'^3T'') 


<»;>    ,   <2i-\> 


1    ^,)i-^/      <;?;>    .   <iii>\  v~^       /    <'>j-l>r    ■   <'">    .  <li-l>        ,   <in>    ,   <sj-i>-< 

^ll    "'    i^Xa      ^3,/J      )+     E      '^  t^3,<.       'l'3,/3  +^3,/3      ^3,a  ] 

1=1 

J 

(m-l)/2 


2Yc/j     =-['\'xa    "^xe    J 


oY.T'  =|[C-"3":r-/' ^'^^c -/'"'"■-"" ^ ^  ^  /,<2-i>,.<^-i> 


C"]-  E  /'^"■""^3T 


r=l 


<;»>         1    ,  <in> 
lYa3        =7t3,a 

0Y33      =  ^3  (2.56b) 

Expressions  (2.56b)  define  quantities  for  layers  above  the  middle  plane  (superscript 
m  is  odd).  The  quantities  for  layers  below  the  middle  plane  (superscript  m  is  even) 
are  given  by  expressions  (2.56c)  which  are  obtained  from  (2.56b)  by  implementing 
the  following:  (1)  replace  li"""  with  -h"^""  where  /z>0;  (2)  change  the  limits  on  the 
summations  to  {r  =  2  to  m/l);  and  (3)  replace  the  summation  index  2/--1  by  2r-2. 
From  (2.55)  and  the  stress-strain  relations  (2.28),  we  obtain 


•^a/J      "  '-aPyn  lolyn        ^  lYyM        ^  2Yy/x     J       *-'a/333   0Y33 

5<;>        ^^<r>   r       <r>         <,•>         <'>  i 
a3      =2Q3/33LY/J3     +^  iY^3     1 


<r>         <r> 


'33 


33a^  LYa/3     +^  lYa^     +^  jYa/J    ]  +  ^333    0Y33 


(2.57) 
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Finally,  from  (2.57)  and  (2.51),  we  obtain  the  following  constitutive  relations 

■'^'a/3      ~l^a/3YM    o^ym      *T*OiPyn    lYy/i      *t^aPyti   2^ yji     "^l-°a/333   0Y33 


-t^a/?       ~2"aPyii    oYyM      "^S^o/JyM    I^yM      *4^aPyli    Z^yli      "*'2^a/333    0^33 


,  <)•>  0  >  r.  <'■>  <'  > 


a         =  2  U^a3y33   oY/33     +l^a3/33   lY/33     J 

rp<r>        ^  r    r)</>  <r>        n  <'>  '^''^  1 

^a         -  "  L  A3;33   oYj33     -^2^03/33  lY;33     J 

^33       -o-°33a/3    oYa^     +1^33a^    lYa^     +2^33a/3   2Ya^     •*'o^3333    0Y33 


where 


+/! 


nr^a 


0yn 


Qym^"'  rfz^^^  ,  (/n  =0,1,2,3,4) 


-/)"'' 

+/!*'^ 


5 


a/333 


^  C:;i,z<^^"'dz<'-^  ,  (m=0,l,2) 


-/i"" 
^h"'^ 


B 


a3Y3 


rj<r> 

/^'"ssym 


jy<r> 

0-^3333 


-/r 


C3?3^"'"'^^"'  >  (m=0,l,2) 


C3i;;2-->"'dz<'-^  ,(m=0,l,2) 


*-3333  "^ 


(2.58) 


(2.59) 


Since  for  the  isothermal  case,  the  elements  Q^  are  constants,  thus  all  ,„5yx/  ^'^'th  odd 
m  will  vanish. 


CHAPTER  III 
THEORY  OF  THREE-LAYERED  PLATES  AND  SPECIAL  CASES 


3.1  Variational  Integral  for  the  Three-Layered  Plate 

The  application  of  Hamilton's  principle  to  thick,  homogeneous  and  isotropic 
plates  using  nonlinear  strain-displacement  relations  has,  according  to  Ebcioglu  [21], 
been  shown  by  Herrmann  and  Armenakas.  It  has  been  used  by  Yu  [15],  for  the  case 
of  small  deformation  and  small  angles  of  rotation,  to  derive  the  equations  of  motion 
of  sandwich  plates.  It  has  been  also  shown  in  [21]  to  lead  to  the  equations  of  motion 
and  boundary  conditions  of  sandwich  plates  for  the  case  of  large  displacements  and 
large  angles  of  rotation.  Ensuingly,  the  variational  integral  for  a  three-layered  plate 
is  derived  from  Hamilton's  principle. 

The  volume  integral  of  (2.48)  is  divided  into  three  parts  corresponding  to  the 
three  layers  of  the  plate.  The  surface  integral  is  divided  into  five  parts  corresponding 
to  the  upper  face,  lower  face,  and  the  edge  surface  of  each  layer.  In  each  of  the 
eight  integrals  the  nine  displacement  components  are  substituted  using  (2.35),  and 
are  now  expressed  in  terms  of  the  twelve  displacement  functions  (2.36).  With  the  use 
of  (2.50),  an  integration  in  the  transverse  coordinate  through  the  thickness  of  each 
layer  is  possible  in  the  three  volume  integrals  as  well  as  in  the  three  surface  integrals 
corresponding  to  the  edge  surface  of  each  layer.  The  result  is  five  surface  integrals 
and  three  line  integrals.  The  surface  integrals,  thus,  are  evaluated  at  the  middle 
plane  of  the  respective  layer.   However,  based  on  the  first  set  of  relations  in  (2.29), 
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the  surface  integrals  are  grouped  and  evaluated  over  the  middle  plane  of  the  middle 

layer,  A^.    Similarly,  the  three  line  integrals,  based  on  the  first  set  of  relations  in 

(2.29)  and  since  Ji^^  =0^^,  are  grouped  and  evaluated  around  the  middle  plane  of 

the  middle  layer.  It  should  be  emphasized  that  the  stress  continuity  conditions  at  the 

interfaces  of  the  layers  are  completely  satisfied  in  the  variational  integral.     In 

addition,  the  transverse  shear  stress,  S^,  at  the  upper  and  lower  faces  vanishes 

automatically,  hence  satisfying  the  condition  for  parabolic  distribution  of  the  shear 

stress,  which  vanishes  on  the   surfaces.     This  condition  is  usually  enforced  by 

additional  equations  in  the  derivation  of  the  variational  integral  of  a  layered  plate 

when  a  cubic  displacement  distribution  in  the  thickness  coordinate  is  assumed,  as  in 

Reddy's  derivation  [29]. 

Define   the   following   notation   for   simplicity   of  terms   recognition,   the 

superscript "/",  "o"  or  "//"  associated  with  a  field  quantity  identifies  its  affiliation  to  the 

upper  layer,  middle  layer  (core)  and  lower  layer,  respectively.    Therefore,  upon 

implementing  a  lengthy  mathematical  derivation,  carrying  out  the  integrations 

outlined  in  the  preceding  paragraph,  and  grouping  of  terms,  the  variational  integral 

for  the  three-layered  plate  is 
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(3.1) 


Using  the  fundamental  lemma  of  the  calculus  of  variations,  we  obtain  twelve 
equations  of  motion  and  appropriate  boundary  conditions  for  any  given  shape  of 
three-layered  plate. 

The  stress,  body  force,  acceleration  and  edge  surface  force  resultants  defined 
in  section  2.6  hold  valid  for  the  present  case  when  the  appropriate  superscript  is 
used.  The  stress  resultants  appearing  in  (3.1)  are  given  by  the  constitutive  relations 
(2.58)  with  the  appropriate  superscript.   Expressions  (2.56)  reduce  to 

o  I  t      o  o  0        0^ 

lYa/)     =  Y  [  ^a,p  +  4f^.a  +"3,«  l|'3,/3  +"3,^  'I's.a  ] 
o  1  r    ,  o       ,  0     ■, 

2YC./J  =-_[^3,a'\'xp]  I  (3.2a) 

O  I    r     ,   O  On 

0Yct3    =7['l'a+"3,a] 

o  1,0 

lYa3    =7^3,a 

o  ,  o 

0Y33    =  ^3 


oYa^   =  ^[{lia   ^It"  K+li  '  Va)^p  ^(th  ^fl"  Vp  ^h'  ^ll'p),^ 
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"  1  r    1  "       1  "     1 

0Y:3  4[^M"/-/^"^3-/^"t3)J 

lYa3=  T'l's.a 

ii  ,  ii 

0Y33  =  ^3 

and  must  be  substituted  into  (2.58)  to  obtain  tlie  stress  resultants. 


(3.2c) 


3.2  KirchhofFs  Hypothesis 

Approximate  theory,  also  known  as  Kirchhoff-Love  plate  rests  on  fundamental 
assumptions  which  are  direct  results  of  Kirchhoffs  Hypotheses.  These  were  first 
introduced  by  Kirchhoff  and  later  presented  by  Love  and  extended  to  thin  shells. 
The  Kirchhoff-Love  plate  is  characterized  by  the  following  assumptions:  (1) 
deflection  of  mid-plane  is  small  compared  to  thickness;  (2)  mid-plane  remains 
unstrained;  hence  mid-plane  does  not  experience  any  in-plane  displacements;  (3) 
plane  sections  normal  to  mid-plane  remain  plane  after  deformation;  thus  the  in-plane 
displacements  are  assumed  to  be  linear  functions  of  the  transverse  coordinate;  (4) 
plane  sections  initially  normal  to  mid-plane  remain  plane  and  normal  to  that  surface 
after  bending;  thus  transverse  shear  strains  are  negligible  and  consequently  deflection 
of  the  plate  is  associated  principally  with  bending  strains,  and,  therefore,  it  is  also 
assumed  that  the  normal  strain  is  negligible;  and  (5)  normal  stress  to  mid-plane  is 
small  compared  tt)  other  stress  components  and  may  be  neglected.  It  is  argued  thai 
the  fifth  assumption  produces  an  inconsistency  in  the  theory.  This  inconsistency  will 
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become  evident  in  the  following  section  in  deriving  von  Karman  equations  and  will 

be  examined  then.   However,  it  is  appropriate  to  clarify  this  point  by  presenting  the 

following  argument  by  Boresi  and  Sidebottom  [41,  pp  452],  which  best  describes  this 

dilemma: 

In  conventional  plate  theory,  it  is  assumed  that  the  plate  is  in  a  state  of  plane 
stress;  that  is,  o^  =  o^  =  a^  =  0.  For  isotropic  elastic  planes,  the  relations 
o^  =  a  =0  are  consistent  with  Kirchhoff  approximation,  which  signifies  that 
e^=€  =0.  However,  the  Kirchhoff  approximation  has  been  criticized  since 
it  includes  the  approximation  e^  =  0.  The  condition  e^  =  0  conflicts  with  the 
assumption  that  a^  =  0.  The  condition  e^^^O  is  incorrect;  however,  the  strain 
e^  has  little  effect  on  the  strains  e^^,  e^,  e^^.  Thus,  the  approximation  £^^  =  0 
is  merely  expedient.  In  the  stress-strain  relations,  the  condition  of  plane  stress 
CT^  =  0  is  commonly  used  instead  of  e^  =  0,  and  this  circumstance  is  often 
regarded  as  an  inconsistency.  However,  in  approximations,  the  significant 
question  is  not  the  consistency  of  the  assumptions,  but  rather  the  magnitude 
of  the  error  that  results,  since  nearly  all  approximations  lead  to 
inconsistencies.  In  plate  theory,  the  values  of  e^^  and  o^^  are  not  of  particular 
importance.  Viewed  in  this  light,  the  Kirchhoff  approximation  merely  implies 
that  e^^  has  small  effects  upon  o^^  and  o^,  and  that  a^^  and  a^^  are  not  very 
significant.  We  observe  further  that  the  Kirchhoff  approximation  need  not  be 
restricted  to  linearly  elastic  plates;  it  is  also  applicable  to  studies  of  plasticity 
and  creep  of  plates,  and  it  is  not  restricted  to  small  displacements. 

Let  the  two  face  layers,  upper  and  lower  layers,  be  characterized  by  only  the 

fourth  assumption  of  those  stated  in  the  previous  paragraph.  Hence,  the  assumptions 

Ya3  =  0  and  Y33~0  when  applied  to  the  displacement  components  (2.35)  and  noting 

the  displacement  continuity  conditions,  leads  to  the  relations 


Vs^O,         11/3^0 


(3.3) 


When  conditions  (3.3)  are  applied  to  the  variational  integral  (3.1),  and  use  is 
made  of  the  Generalized  Gauss's  theorem  and  integration  by  parts,  i.e., 
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St.SVnds 


Si,^  SVdA 


(3.4) 


where  Si  is  a  second  order  tensor,  we  obtain 
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3.3  Von  Karman's  Equations  of  Plate  Large  Deflection 

The  equations  of  equilibrium  and  boundary  conditions  for  the  large  deflection 
of  isotropic  elastic  plates  are  derived  from  the  variational  integrals  (3.1)  by 
eliminating  the  core  and  assuming  identical  faces  (geometrically  and  materially), 

fi"  ^  0   and   2h'  =  2h"  =  t/2  (3.6) 

Also,  the  following  assumptions  are  imposed:  (1)  Kirchhoff  hypothesis  holds;  (2) 
transverse  deflection  is  uniform  through  the  thickness  of  the  plate,  hence  i|fj— ^  0  and 
i|r"  — t  0;  and  (3)  static  case  and  zero  body  forces.  Therefore,  from  (2.35)  and  (3.3), 
the  deformation  of  the  plate  is  now  described  by 


jrii  o       /     a      ,   ,;  >.       o  jrii  o 

Hence,  V^  =  u°  - z u^,^       ,      V^  =  u^ 


-li'<z'<  +li' 
-h"<z"<  +h'' 
■2h"<  z  <  +2//' 


(3.7) 


where,  w°  and  w^  describe  the  in-plane  and  transverse  displacements  of  the  middle 
plane  (z  =  0),  respectively.  We  also  assume  the  following  loading  conditions  on  the 
upper  and  lower  faces, 


f 


t'  =  0    ,       t"  =  p 

0~-k  '  0-1.  -r 


(3.8) 


Upon  imposing  these  conditions  into  the  general  variational  integral  (3.1),  we  obtain 


(3.9) 
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The  definitions  (3.2)  reduce  to 
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And  the  nonvanishing  coefficients  in  (2.59)  are 
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(3.11) 


Hence,  in  view  of  (3.10)  and  (3.11),  and  from  the  constitutive  relations  (2.58),  it  can 
be  easily  shown  that 
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(3.12) 


In  arriving  at  the  correct  form  for  the  two  dimensional  fourth  order  isotropic  tensor 
^aPyiL^  we  must  derive  it  from  the  complete  three  dimensional  tensor  by  imposing 
the  conditions  of  the  Kirchhoff  hypothesis  on  the  strains.  However,  this  process  will 
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lead  to  the  plane  strain  case,  and  consequently  the  flexural  rigidity  constant  D  will 
differ  from  that  of  von  Karman,  since  he  assumed  the  plane  stress  case.  The 
difference  is  demonstrated  by  the  ratio 
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(3.13) 


For  the  plane  stress  case, 
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We  substitute  (3.12)  and  (3.14)  into  (3.9)  and  use  the  fundamental  lemma  of  the 
calculus  of  variations  to  obtain  the  following  equations  of  equilibrium  and  boundary 
conditions  for  the  von  Karman  plate  type. 
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D  = 


P  E 


12  (1-U-) 


and 


E  =  2G(l+u) 


(3.16) 


and  E,  G  and  u  are  the  elastic  modulus,  shear  modulus  and  Poisson's  ratio, 
respectively.  The  stress  resultants  N^p  are  due  to  the  middle  plane  in-plane  strains, 
and  they  are  unknown  quantities.  A  stress  function,  0  =  cp{x^)  is  introduced,  such 
that  the  first  two  equilibrium  equations  in  (3.15)  are  satisfied  identically  by  0(.Va), 
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and  the  stress  function  is  related  to  the  resultants  as  follows; 

N,p  =  t{6,p<p,^^-4>,^p)  (3.17) 

Two  equations  of  equilibrium  are  thus  obtained  in  terms  of  the  stress  function  (pix^) 
and  the  vertical  deflection  of  the  middle  plane,  u^;  and  they  are  identical  to  those  of 
von  Karman.    When  written  in  expanded  form,  these  are 

"^^      dx-d^'-  ^~^      D     1*^^      dx'~  dy-         dxd)>  dxd}' 

d^cp       <?*0  _  9      d^(p      -  f  \  (  ^^  ^2  _  ^^v  d^^  -I 
dr^     ^     ^  dx-dy"  dxdy         dx'   dy- 

The  first  equation  in  (3.18)  is  obtained  from  the  third  equilibrium  equation  in  (3.15) 
by  substituting  (3.17).  The  second  equation  in  (3.18)  is  obtained  by  enforcing  the 
strain  compatibility  condition  which  is  violated  by  the  assumptions  on  the  transverse 
strain  components. 


CHAPTER  IV 

THE  SOLUTION  TO  THE  SANDWICH  PLATE 

BENDING  AND  BUCKLING  EQUATIONS 


4.1  The  Assumptions  of  Eringen's  Equations 

The  problem  of  bending  and  buckling  of  sandwich  plates  has  been  formulated 
and  solved  by  Eringen  [11].  The  underlying  assumptions  for  the  problem  formulation 
are  summarized  as  follows:  (1)  usual  small  deformation  theory;  (2)  faces  are  very  thin 
as  compared  to  the  core,  but  are  not  membrane;  (3)  middle  plane  of  the  plate 
remains  unstretched  subsequent  to  loading;  (4)  displacements  of  the  plate  are  linear 
functions  of  the  distance  from  the  middle  plane  of  the  plate;  and  (5)  linear  theory 
of  elasticity.  The  problem  formulation  considers  all  sbc  components  of  the  stress 
tensor  in  the  core  and  the  bending  rigidity  of  the  face  layers. 

The  simplicity  of  Eringen's  theory  is  due  mainly  to  its  linearity  while  some  of 
the  generality  of  the  problem  is  retained.  The  assumption  of  a  linear  distribution  of 
the  deflection  components  does  not  extend  through  the  entire  plate  thickness.  More 
specifically,  it  is  defined  only  up  to  the  middle  planes  of  the  face  layers.  This 
consideration  limits  the  range  of  application  of  Eringen's  theory  to  sandwich  plates 
of  very  thin  faces  as  compared  to  the  core. 

4.2  The  Equations  of  Equilibrium  and  Boundary  Conditions 

The  general  variational  integral  and  all  fundamental  equations  for  the  three- 
layered  plate  were  derived  in  the  previous  chapter.  For  the  static  case,  with  no  body 
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forces  present,  the  equations  of  equilibrium  and  the  necessary  and  sufficient 
boundary  conditions  are  obtained  from  the  variational  integral  (3.1)  using  the 
fundamental  lemma  of  the  calculus  of  variations.  In  order  to  make  a  comprehensive 
and  relevant  comparison  of  the  results,  the  face  layers  are  assumed  to  be  identical; 
hence  h'  =  If.  Kirchhoffs  hypothesis  is  assumed  valid  in  the  face  layers;  thus 
conditions  (3.3)  are  valid.  In  addition,  we  assume  that  the  middle  plane  of  the  plate 
(which  is  also  that  of  the  core  layer)  is  unstretched  subsequent  to  deformation,  hence 
ul  =  0.  The  above  mentioned  assumptions  allow  for  the  comparison  with  Eringen's 
results;  however,  while  Eringen's  problem  is  based  on  a  two-dimensional  linear 
elasticity  formulation,  the  present  is  based  on  a  three-dimensional  non-linear 
elasticity  formulation.  The  major  and  underlying  difference  is  that  the  present 
formulation  assumes  the  displacement  distribution  to  be  linear  in  each  layer  rather 
than  through  the  entire  plate.  Implementing  the  aforementioned  assumptions  into 
(2.35),  we  obtain 


(4.1) 


which  describe  the  displacement  distribution  in  the  three-layered  plate. 

The  loading  on  the  upper  and  lower  faces  of  the  plate  is  assumed  to  be; 

o^a    =     oC    =     4   =  0     ,       4   =  p{xy)  (4.2) 

From  (3.1),  (3.3),  (4.1)  and  (4.2)  we  obtain 
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M°  K° 

h  h° 

rpO  KT^ 
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.  { (S;  .5;  .5;')  -  [  (<  .<,  .<,  )<„  .  h  ^N^^  .  ^  -yvj; )  ^l^ 

h 

/z  h  ii " 


^  {Ml,  ^h  -Nl,  -M:,  +/z  'WJi  ),^  +  Zi  ]^«p  /z "  5 Tlr"  )  ^.  =  0  (4.3) 


The  assumption  u°^  =  0,  according  to  Ebcioglu  [22],  allows  us  to  linearize  the 
strain-displacement  relations.  Therefore,  we  can  omit  N^,,  Ml,  and  Kl,  in  all 
non-linear  terms.  Also,  in  view  of  (3.3),  (2.57)  and  (2.59),  the  necessary  stress 
resultants  appearing  in  (4.3)  are  obtained  from  (2.58): 
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M<;,  =  -4/i'  C;,,,[a3%/i°il;^],^, 


/l" 


<  =  -/^"/^"c;,,,  [^l,-K,,-2!j-^iu,"-h'^  ¥,),,,] 


K  =  -^h-''c:p^^[u^-h'^^]fl],^^ 


where,  Q^^    =  Q^^  ,  /r'  =  /i"  and  (3.14)  is  valid. 

By  substituting  (4.4)  and  (3.14)  into  (4.3),  we  obtain 


(4.4) 


/^  O  T  C 

4/;'G'{[l  +  ^ 

A.    \  Qiji 


Vi^^sp-l 


i+v'        cr 


2h' 


■]h''re,,p--^^^—u 


XaPP 


1-v'     G7'(1-2v'')"         "'"^      (1-v') 

^/.>:-l^<}/.«5C  .  {[(A^:,-/vp<.(yv;^-A^;)/.°<j,, 

1  Lr  1 


^E'h'  ,0,0 


1-v'^ 


h°^lpaa-^D,ul,,^^  ^ICh^^ll,  ^ICh^ulpp  -p)Su° 


^  ^^Kp  -Kp)<^^.p  -C,)/^  >3,a],,  -  8D,/^>3,.a/J^ 


.i£!^(/r<.„-^^D-p)/^^^t3V"  ^i  v^^:^^-^^^ 


/' 


9  /^  °  /)  " 


L° 


8G7i'v'     4G''/z%'' 

+ 


1-v'        3(l-2v'')J 
1-v' 


t^'^.Kp 
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{(L;  -Ll^-h'iS^  -Sp  -  [4G'Wz°(il;:^  -ij;;.,  --^C5„^) 


1-v' 
1  -v' 

-(i!^;   -Nf^)/i>°,  -8D^/i>^,,,^  -^^I^il;°^;5],n^}/z^         )  ds  =  0         (4.5) 


where; 


^1  = 


2£'7z' 


^„^20i^^        jo^2hy^        j.^^j^^l^o^    (4.6) 
3(l-v'^)  (l-2v°)  3 

By  using  the  fundamental  lemma  of  the  calculus  of  variations,  we  obtain  from 

(4.5)  the  following  four  partial  differential  equations  of  equilibrium  in  cartesian 

coordinates; 

8£>i<aaw-2/z''G°«,,-iI;;,^)-i^/z''<,^^ 

(1-v') 


8Dj/r<'H;°,,„^^ 


4G°I°h' 


r 


I  o  •  o         J  a     ,  o 

1^     ^3,a«  -  — -  ^3 

h 


=  -P  -  {^0  -  Kp )  «^.^  -  (N^p  -  < )  li "  ^l.p 


r^  o  1  o 


l  +  v'      G"r{a°-\) 


1-v' 


G'/' 


^"^^30. 


lG°h° 


2h' 


(  4fa  +  «3,a  )  -  ■  ^3aP0  =  ^ 


G'/'  1-v' 

and  the  boundary  conditions  are  specified  from 


(4.7) 
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5il;°=0,      or      5'  -5"  +  -^ 

-a       -a        f^o 


(4G'h''^^^^){h'>Va,p^ll''Vp,J 


8G'h'v'      ACh^v" 

+ 


[      1-v'  3(l-2v'')J 


K,  h "  V,, 


-8G'7i'"«„^-_^5„^<,,) 


1-v 


o^P 


1-v 


SuL=Q,      or     L'  +L"+h'iS'  -S")  = 

-8Dj(l-v')«„^  +— ^^^/j^s^yy) 

1-v' 


o'H 


5il;°    =0,      or     L'- 


1-v' 


5^3-0,      or     5'+5°+S" 
■*  —3—3     —3 


(A^,  -A^)"i  MA^,  -K.Vf^l,  -SD.ula, 


'  i  U  ii'  U  o  /  ,,0  1  +  V'   _,_o 


+4G'/z"/z°(Tl;;,,, +-L_r-il;:.„^)-^2G''/r(t^+f/;^) 

1-v' 


o''^ 


<Silf3=0,      or     S'+-J--5"  = 
-3     ,^.     -3 


-8Z)i/^"il/U,-^^^t3,, 


(4.8) 


.'^/! 


The  above  partial  differential  equations  are  comparable  to  those  derived  by 
Ebcioglu  [22]  under  the  same  assumptions.  They  represent  four  coupled  but  linear 
fourth  order  partial  differential  equations  to  be  solved  for  the  four  unknown 
displacement  functions 

("3^^3, 4^:)  (4.9) 

The  boundary  conditions  are  given  in  (4.8)  and  must  be  satisfied  by  the  solution. 
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Three  special  cases  of  the  above  equations  are  derived  by  Ebcioglu  [22]  for: 

(a)  no  core,  h°  -►O; 

(b)  bending  moment  of  the  core  is  negHgible,  D°  -  0;  and 

(c)  D°  -►O  and  i|f^  -♦O,  which  leads  to  the  equations  given  by  Hoff  [9]. 

4.3  Applications:  Numerical  Results  and  Discussion 

Consider  a  rectangular  sandwich  plate  simply  supported  along  face  edges  x  =  0, 
x=a  andy  =  0,y=6.  It  is  submitted  to  an  arbitrarily  distributed  transverse  load/7(x^) 
on  its  face  and  a  uniformly  distributed  axial  compressive  load  {-N^)  along  the  face 
edges  x  =  0  and  x=a.  The  problem  is  solved  by  solving  the  four  partial  differential 
equations  (4.7)  with  N'^=N.l  =  -N,  and  A^^'^  =  A^^'^  =  0,.  The  boundary 
conditions  specified  in  (4.8)  become 

along  x  =  0  andx=a: 

O  ,  O  ,  O  /-, 

G  h  oG  h         ,  .0       i  ,  o  \         V       /  I  o    ,  ,  o  s.      ,  o 

G"h"{l-v')        '  '        l-2v'' 


G'h' 


("3,11 +  v  "3,22)  =0 

1-v' 


—  (^fi.i  +  V  11/2,2)  -^i  ("3,11  +  V  1/322)  =  0 

1-v' 

along  y  =  0  and  y=b: 
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O  ,0  ,  O         f^ 


GVi 


6G'li 


G"h°il-v') 


(  4^2  2  -  V'  ^u  )  ^  — ^  (  ^1,1  "  <2)  -  t2,2 

l-2v° 

G    h       /      o  i      o     \       f\ 

-  -(«3,22  +  v'i(3ii)  =0 

1-v' 


G'fi'  h' 


(  ^2,2  +  V  ti,!  )  -  Dj  (U^22  +  V'  Wa*;!!  )  =  0 


1-v' 


(4.10b) 


It  may  be  seen  that  the  following  series  satisfy  all  of  the  boundary  conditions: 


CO  00 


u^{x,y)  =  Y.   E  "^' 


3;?!/! 


sin 


mirx     .    niry 


sin- 


m=\     n=\ 


O)  00 


a 


'^°i{x,y)  =  Y,   E   ^3«m  sin 


miTX    .   niry 
sin 


m  =  l     n-l 


00  CO 


ti(j^,>')  =  E  E  ^1'""  ^os 


m  =  \     n  =  \ 


CO  00 


miTx     ■    mry 
sin — - 

a  b 


mirx         niry 


(4.11) 


^2i^,y)  =  Y.     E     ^2.„«   sin COS 

where  u^„„,  i|;3„„,,  il/^,,,,  and  i|j^„„  are  unknown  constants. 

The  foregoing  series  satisfy  differential  equations  (4.7)  if  p(;c^)  is  expanded 
into  the  Fourier  series 

mitx  ^;„nny  (4.12) 


P(x,y)  =  Yl     E    Pmn    Sin- 


sin- 


m=l     71=1 


a 


where 


P     =- 

ab  J 


0 


/      .    .mirx     .    niry    ,    , 
p{x,y)  sm sin— ;-^  ax  ay 


(4.13) 


a 


Substitution  of  (4.11)  and  (4.12)  into  equations  (4.7)  yields  four  simultaneous 
algebraic  equations  for  the  unknown  coefficients  u°,„„,  ^l,„„,  Vh,,,,  ^^^  ^Iw,  (see 
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Appendix  B).   The  displacement  functions  for  any  point  (x^)  can  thus  be  obtained 

by  means  of  series  (4.11),  substituted  into  (4.1)  to  obtain  the  components  of 

displacement  at  any  point  of  the  3-D  space  of  the  plate. 

The  equations  presented  by  Eringen  are  solved  for  the  deflections  of  the 

sandwich  plate  in  the  case  of  a  simply  supported  rectangular  plate  under  the  actions 

of  a  uniformly  distributed  transverse  load  on  the  upper  face  of  the  sandwich  plate 

and  a  uniformly  distributed  edge  compressive  load  along  the  edges  x  =  0  and  x=a. 

In  this  cSiSQ  p(x^)  =  p  =  constant  and  equation  (4.13)  reduces  to 

P_  =  J^,  m,Az  =  1,3,5,..  (4.14) 

IT   m  n 


The  deflections  of  the  middle  plane  are  computed  based  on  two  approaches, 
Eringen's  theory  and  the  present,  for  several  different  geometrical  and  material 
configurations.  While  maintaining  as  constants  the  following:  (p  =  100  psi;  N^  =  0; 
a  =  100  inches;  £^=2*10'  ksi;  and  h  =t,  +  2tf=  10  inches,  where  t,=2/i°  and  tf=2h')  and 
with  b/a  - 1,  several  analyses  are  conducted  by  varying  the  ratios  E^/E^  and  t^/t^.  The 
number  of  Fourier  terms  in  the  series  was  taken  to  be  19.  The  results  of  these 
analyses  are  conveniently  plotted  in  comparable  groups  in  Figures  4.1  through  4.4. 
The  deflections  of  the  middle  plane  of  the  plate  are  plotted  for  half  the  span  length 
in  the  x-direction  along  the  line  y= 6/2,  which,  due  to  the  loading  and  geometrical 
symmetry  of  the  problem,  is  representative  of  the  deflected  shape  of  middle  plane 
of  the  plate  in  other  perpendicular  directions.  In  Figure  4.1,  the  face  to  core 
thickness  ratio  was  held  at  1  to  10  which  indicates  very  thin  faces  as  compared  to  the 
core.  The  deflected  shape  predicted  by  the  two  approaches  for  an  elastic  modulus 
ratio  of  1  to  10  are  almost  identical;  however,  as  the  elastic  modulus  ratio  increases, 
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Distance  in  x-direction  along  line  y=b/2 
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indicating  a  weaker  core  material,  Eringen's  results  are  slightly  higher  than  the 
present  results.  The  maximum  difference  at  the  center  of  the  plate  is  of  the  order 
of  9%  which  is  still  within  an  acceptable  comparison  range.  Figures  4.2  and  4.3 
present  the  results  for  face  to  core  thickness  ratios  of  1  to  4  and  1  to  2,  respectively. 
The  pattern  is  similar  to  that  in  Figure  4.1,  where  again  the  difference  between  the 
two  predicted  deflected  shapes  tends  to  increase  as  the  core  becomes  weaker.  The 
variation  between  the  two  predicted  deflected  shapes  seems  to  widen  with  the 
increase  of  the  face  to  core  thickness  ratio  as  well.  The  results  become  incomparable 
when  the  thickness  ratio  becomes  unity,  that  is  the  face  and  core  layers  are  of  equal 
thicknesses. 

The  results  of  Eringen's  equations  start  diverging  from  the  present  results  and 
become  unreliable  by  over-predicting  the  deflected  shape  enormously  as  the  face 
thickness  approaches  that  of  the  core,  as  well  as  when  the  core  material  tends  to 
weaken.  In  order  to  justify  the  just  stated  conclusion,  consider  the  following 
situation.  The  three-layered  plate  in  the  previous  analysis  is  subjected  to  the 
following  {b/a  =  10,  a/h  =  1,  and Ef/E^=\)  which  simply  states  that  the  face  and  core 
layers  are  assumed  to  be  made  of  the  same  material.  By  holding  the  thickness  of  the 
plate  constant,  the  results  may  be  compared  to  that  predicted  by  the  classical  plate 
theory.  The  maximum  deflections  at  the  center  of  the  plate  are  computed  based  on 
the  present  formulation  and  based  on  Eringen's  equations  for  several  ratios  of  face 
to  total  thickness  of  the  plate  and  the  results  are  plotted  in  Figure  4.5  for  the 
maximum  deflection  versus  the  thickness  ratio.  The  values  compare  very  well  for 
small  thickness  ratios;  however  as  the  ratio  increases,  more  specifically  beyond  1/3, 
the  maximum  deflection  predicted  by  Eringen's  equations  increases  rapidly.  The 
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Figure  4.5:  Comparison  of  vertical  deflection  at  the  center  of  a  simply  supported 
square  plate  subject  to  a  uniformly-distributed  transverse  load. 


maximum  deflections  for  a  homogenous  plate  of  the  same  dimensions  and  thickness 
are  obtained  from  the  equations  of  the  classical  plate  theory,  and  the  ratios  of  the 
maximum  deflections  from  the  present  and  from  Eringen's  equations  to  that  of  the 
classical  plate  theory  are  plotted  versus  the  thickness  ratio  in  Figure  4.6.  The  ratios 
of  the  present  formulation  to  the  CPT  deflections  are  very  close  to  unity,  indicating 
a  good  prediction  by  the  present  equations  regardless  of  face  to  total  thickness  ratio. 
They  are  not  exactly  equal  to  unity  for  the  simple  reason  that  in  the  present 
formulation  the  transverse  shear  strains  are  not  neglected.  The  ratio  of  Eringen's  to 
CPT  deflection  starts  diverging  from  unity  when  the  face  to  total  thickness  ratio 
reaches  1/3. 
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(Note:  MLPT  =  Multi-Layered  Plate  Theory;  CPT  =  Classical  Plate  Theory.) 


4.4  Stability  Analysis  of  the  Sandwich  Plate 

The  transverse  displacement  of  the  plate  is  assumed  to  vary  linearly  in  the  z- 
coordinates,  thereby  allowing  the  core  to  experience  a  change  in  thickness.  This  is 
usually  referred  to  as  the  flattening  of  the  core  (i.e.,  the  approach  of  the  two  faces). 
It  is  measured  as  the  difference  between  the  transverse  displacements  of  the  two  face 
layers.  It  was  first  examined  by  Eringen  [11]  as  to  its  effect  on  the  overall  instability 
of  sandwich  plates.   The  present  theory  predicts  the  flattening  of  the  core,  and  the 
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results  are  in  good  agreement  with  Eringen's  for  the  same  ranges  stated  earlier  in  the 

case  of  overall  plate  deflections. 

For  the  particular  example  of  an  arbitrarily  distributed  transverse  load,  it  can 
be  easily  proved  that,  owing  to  the  in-plane  compressive  forces  A''_^,  all  the  coefficients 
w|„,„,  as  well  as  all  the  coefficients  i|;|„„„  increase;  hence  deflections  of  a  compressed 
plate  are  larger  than  those  of  the  equally  loaded  identical  plate  without  in-plane 
compressive  loads  applied  at  the  face  edges.  It  is  seen  also  that  by  a  gradual 
increase  of  compression  we  approach  a  value  of  N^  for  which  one  of  the  coefficients 
u°„,„  or  one  of  the  coefficients  i|rj'„,„  approaches  infinity.  The  smallest  of  these  values 
of  N^  is  called  the  critical  value. 

From  the  solution  to  the  four  simultaneous  algebraic  equations  for  the 

unknown  coefficients  u^„„,  ^^l,„„  ^]f1„„  and  ij;^,,,,,  we  obtain  expressions  for  the  two 

coefficients  w|,„„  and  i|f|„,„  as  the  following  (see  Appendix  B) 

P  P 

o  '■  inn  ,  o  mn  (A  1  <\ 

mn  mn 

Here,  buckling  occurs  when  either  yl^„-+0  or  5^„-0,  depending  upon  whichever  gives 
the  smallest  critical  load.  These  critical  values  are  N^nj  and  A^c;?.2'  ^"^  ^^ey  are 
positive  increasing  functions  of  n.  Thus,  n  =  \  makes  A^qj.^  and  Ncr,2  minima.  The 
critical  stresses  a^/^;  and  CTc^2  may  be  obtained  by  dividing  A^^i?.;  ^i^^  ^cm  by  face 
thickness  (2/z').  The  critical  stress  acR,i  represents  the  buckling  of  the  plate  as  a 
whole  and  o^ni  represents  the  buckling  due  to  flattening  of  the  core. 

To  better  demonstrate  the  effect  of  the  compressibility  of  the  core  (see  Figure 
4.7),  consider  the  sandwich  plate  having  faces  made  of  aluminum  alloy  materials 
(£y=  10'' psi)  and  the  core  made  of  soft  isotropic  materials   (£'^  =  2000  psi).  We 
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EULER  (OVERALL)  BUCKLING 


BUCKLING  WITH  CORE  FLATTENING 


Figure  4.7:  Buckling  modes  of  a  sandwich  plate. 


observe  from  the  computation  of  Ocuj  and  Ocji2  that  (a)  Ocrj  and  Ocrj  have 
minimums  at  some  a/mb,  and  (b)  the  number  of  half-waves  (m)  for  which  (a^/,  ;)n,i„ 
occurs  is  much  smaller  than  that  for  {ocR.2)min- 

The  values  of  (ac^;)^i„  and  iacji2)mm  are  plotted  for  various  face  thickness  to 
width  ratios  tj^/b  as  well  as  for  various  core  thickness  to  width  ratios  tjb,  as  shown 
in  Figure  4.8.  Similar  curves  are  generated  in  Figure  4.9  based  on  Eringen's  theory. 

The  obvious  conclusion  from  Figure  4.8  is  that,  except  for  the  case  of 
t^/b  =  1/1000,  flattening  of  the  core  does  not  influence  the  overall  buckhng  of  the 
plate  until  the  core  thickness  to  width  ratio  exceeds  1  to  5  which  indicates  a  very 
thick  core  of  a  very  narrow  plate  (thus  sandwich  beam-column  analysis  may  be  more 
appropriate).  For  face  thickness  to  width  ratios  of  the  order  of  1  to  1000  and 
smaller,  buckling  occur  due  to  flattening  of  the  core  when  tjb  exceeds  1/10.  These 
results  greatly  under-emphasize  the  effects  of  core  flattening  predicted  by  Eringen 
shown  in  Figure  4.9.  A  statement  to  that  effect  may  be  made,  that  while   the 
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Figure  4.8:  Minimum  Critical  Stresses  of  the  present  theory  for  both  types  of  buckling 
(plate  as  a  whole  &  with  the  flattening  of  the  core). 
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Figure  4.10:  Comparison  of  the  critical  stresses  for  tc/b  =  l/20. 

compressibility  of  the  core  may  be  a  critical  element  in  the  instability  analysis  of 
sandwich  plates,  it  does  not  become  so  until  face  to  core  thickness  ratio  tj-/t^  is  very 
small. 

In  general,  while  predicted  {ocR.2)min  ^^e  tremendously  different,  the  overall 
buckling  or  minimum  o^r  (the  lesser  of  {ocR,i)m\a  ^^^  i°cR.2)mm)  compares  well 
between  the  present  and  Eringen's  theory.  In  Figure  4.10,  these  values  based  on  the 
present  and  Eringen's  theory  are  compared  with  the  minimum  a^y?  from  Reissner's 
theory  [10]  for  core  thickness  to  width  ratio  of  1  to  20.  The  ratio  of  the  values 
increases  as  tJb  increases  or  as  t^/tj-  decreases.   As  indicated  by  Reissner,  this  may 
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be  due  to  the  fact  that  for  thin  sandwich  plates  (small  tjtj)  overall  buckling  occurs 
with  small  wave  lengths  for  which  the  theory  developed  there  becomes  inadequate. 


4.5  Reduction  to  the  Classical  Plate  Bending  and  Buckling  Equations 

When  the  core  is  eliminated  {lf  =  Q),  the  loading  is  assumed  to  be  a  uniformly 
distributed  transverse  load  only,  and,  making  A^^  =  0,  we  may  derive  an  expression  for 
the  vertical  deflections  of  a  homogenous  plate  having  thickness  i  =  4/i'  and  dimensions 
a  and  b.  Implementing  h°  =  Q  into  the  four  simultaneous  algebraic  equations  for  the 
unknown  coefficients  ^3°,,,,  \^l,„„,  i|fi„,„  and  T|;i„„  and  solving  lead  to  the  only  non- 
vanishing  coefficients 
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where 
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TT    m  n 


m,n  =  1,3,5,.. 


The  vertical  deflection  is  obtained  by  substituting  the  coefficients  (4.16)  into 
the  first  of  the  series  (4.11)  and  noting  that  the  second  series  in  (4.11)  vanishes  for 

every  term  in  the  series.  Thus 
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The  deflections  Vs{x^)  given  in  (4.17)  are  identical  with  the  one  obtained  in  the 

classical  plate  theory  [42], 

We  now  assume  that  there  is  a  compressive  force,  A^^.,  acting  on  the  edges  of 
each  of  the  face  layers  (i.e.,  2N^  acting  on  the  plate  edges)  and  that  the  transverse 
load  may  be  arbitrarily  distributed.  Implementing  h°-0  into  the  four  simultaneous 
algebraic  equations  as  before  and  solving  lead  to  the  only  non-vanishing  coefficients; 
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(4.19) 


The  critical  buckling  load  is  the  smallest  value  oiN^  for  which  one  of  the  coefficients 

Wj„„,  becomes  infinite;  thus  it  is  necessary  to  take  n  =  \,  thereby  obtaining 
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which  is  identical  with  the  one  obtained  in  the  classical  theory  [42]. 


CHAPTER  V 
CONCLUDING  REMARKS 


A  variational  theorem  of  three-dimensional  elasticity  has  been  used  in 
developing  the  fundamental  equations  of  the  theory  of  multi-layered  plates  in  terms 
of  a  reference  state.  The  theory  has  been  demonstrated  for  the  case  of  a  three- 
layered  plate  as  far  as  deriving  the  equations  of  motion  and  boundary  conditions. 
Special  cases  were  considered  and  studied  by  solving  the  differential  equations  for 
given  loading  and  geometry  conditions  and  comparing  the  results  obtained 
numerically  to  those  of  the  literature. 

The  theory  presented  considers  all  stress  components  in  each  layer;  and  the 
assumed  displacement  components  were  defined  separately  for  each  layer.  Hence, 
it  can  be  stated  that  the  procedure  of  this  derivation  without  making  any  further 
assumptions  can  be  extended  to  a  layered  plate  of  more  than  three  layers.  The 
mathematical  derivation,  however,  is  more  involved.  Furthermore,  it  may  not  be  a 
worthwhile  task  at  the  present  time  to  attempt  such  a  derivation  unless  means  for 
numerical  implementation  exist.  It  is  obvious  that  the  theoretical  development  in  this 
field  is  more  advanced  than  the  availability  of  techniques  for  solving  the  system  of 
non-linear  partial  differential  equations  characterizing  these  structural  elements. 
This  lack  of  solving  techniques  incorporates  both  analytical  as  well  as  numerical 
aspects.  However,  in  this  age  of  rapid  advancement  in  the  area  of  computational 
utilization,  which  has  yielded  the  introduction   and   implementation   of  numerous 
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sophisticated  numerical  techniques  in  several  areas  of  the  engineering  domain,  it  is 

well  conceivable  to  find  in  the  near  future  some  mechanism  for  solving  these  highly 
non-linear  systems. 

The  general  theory  presented  has  the  advantage  of  producing  a  complete  set 
of  fundamental  equations  consistent  with  various  stages  of  linearization  in  the  general 
strain-displacement  relations.  The  literature  abounds  in  such  intermediate  theories, 
and  the  present  work  is  hoped  to  have  shed  some  light  on  these  special  cases  as  well 
as  cleared  the  way  for  a  systematic  development  of  plate  and  shell  theories  directly 
from  the  three-dimensional  theory  of  elasticity  in  terms  of  a  reference  state. 


APPENDIX  A 
NOTATIONS  OF  SOME  EARLIER  WRITERS 


(a)  Notations  in  general  theory  appearing  in  chapter  II: 


Field  Quantity 

Green  and 

Adkins  [38] 

(Present) 

Fung  [39] 

Malvern 
[40] 

Position  of  undeformed  body 

^i 

«, 

a; 

Position  of  deformed  body 

yi 

-^, 

•V, 

Displacement  components 

K 

", 

", 

Green's  deformation  tensor 

G, 

Q 

Symmetric  Strain  tensor 

y. 

^. 

^. 

Unit  normal  in  undeformed  position 

o'h 

^,n 

Symmetric  stress  tensor  per  unit  area 
of  undeformed  body  (second  Piola- 
Kirchhoff  stress  tensor) 

S'J 

^•J 

T,j 

Unsymmetric  stress  tensor  per  unit 
area  of  undeformed  body  (first  Piola- 
Kirchhoff  stress  tensor) 

fij 

T, 

^  'I 

Stress  vector  per  unit  area  of 
undeformed  body  associated  with 
surface  in  deformed  body 

ot 

t° 

(b)  For  5,  : 


(c)  For  ?,^  : 


Novozhilov 

Eringen 

Yu 

Novozhilov 

Eringen 


a^    {i=x,y,z;j=x,y,z) 

7*' ,  pseudo-stress 

a-  ,  Kirch hoff-Trefftz  stress 

al    {i=x,y,z;  j  =  i,r]X) 

1^  ,  pseudo-stress,  Piola  stress 
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APPENDIX  B 
SOLUTION  TO  PDE  AND  LISTINGS  OF  COMPUTER  PROGRAMS 


Define  the  following  geometric  and  material  constants 


^1  = 


f      \ 

=  6Dj 

h" 

1-v'' 

[h'j 

AG"rh°h''  _  2 

G"h 

o^ 

I' 


V,=h" 


ri,=li° 


1  + 


QOJO 

G'r 


h" 


1  + 


G"h" 
6G'li' 


1+v'  ^G^ria"-!) 


1-v' 


G'l' 


h° 


1+v'     G°h°{a°-\) 


1-v' 


6G'h 


2G"h"  G" 


G'r 

2/2' 


2h'G' 


1-v' 


(B.l) 


The  four  simultaneous  algebraic  equations  for  the  unknown  coefficients  in  Fourier 
series  (4.11)  may  be  written  in  matrix  form  as  follows 

/  \  / 

0        ^22        0  0  U3„j  C, 

1  f     =    " 

$41        0         $43       $44]     [^2n,J  1^4 

where 


(B.2) 


'  1  iitn 


$11  =  8Z), 


C,  =  -F 


^3    =    C4    =    0 


/  miT  v4    ^  .  /HTr  sT  ,  7Z7r  .9     ,  inr  .4 
( )  +2( )'{)-  +  ()' 

a  abb 


+  2/rG" 


a  b 


2NX—? 
a 
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^^3  =  2/j''G°(:^)-r?j 
a 


nil 


b 


?,2  =  ^D^h" 


( )  +(^)(-r-)" 

a  a        b 


,  rrui  ■.I  ,  nil  ..       ,  nn  si 
( )"(-;-)  +  (^) 

abb 


nn  s4 


b  b 
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a  b  i^o^ 
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a  a        b 
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I   III  II     si 
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a        b 


41 
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•44 


ri,i )(— ) 

a         b 


-ri- 


( )"  +  (^)- 

a  b 


nA—-)--r]. 


The  algebraic  system  (B.2)  is  thus  solved  for  the  unknown  constants  (;/-;„,„, 
^3miv  ^imiv  ^2miX  ''^'^'^  thcsc  coustauts  are  in  turn  substituted  into  their  respective 
series  in  (4.11)  to  obtain  an  expression  for  the  displacement  functions  at  any  point 
{x^).  The  constants  (z/,^,,,,,,  i|rj,„„,  i|/;„„„  ^2mn)  ^^e  actually  functions  of  the  series  indices 
m  and  n;  therefore,  the  simultaneous  algebraic  system  of  equations  is  solved  as  many 
times  as  the  number  of  terms  in  the  series.  A  complete  computer  program 
implementing  the  procedure  outlined  in  chapter  IV  and  the  foregoing  paragraph  is 
presented  at  the  end  of  this  section. 

In  order  to  determine  the  critical  load  A^^.  it  is  necessary  to  obtain  an 
expression  for  Wj„,„  and  ^^„,„  and  then  set  the  denominator  equal  to  zero.  It  can  be 
shown  that  this  is  obtained  by  solving  the  following  two  equations  for  N^ 
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i..  =  0 

(^33  ^44  "^34^43) 

where  the  constants  appearing  in  (B.3)  are  defined  above.    Note  that  the  term  A^, 
appears  only  in  the  definitions  of  ijj  and  E,.y  (LHS  of  equations  B.3). 
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Listing  of  a  computer  program  written  in  FORTRAN  language  to  compute  and  print 
the  complete  defomation  data  of  a  simply  supported  rectangular  sandwich  plate 
under  uniformly  distributed  transverse  load,  and  to  evaluate  the  critical  stresses  for 
instability,  based  on  the  present  theory. 

C  Program  to  solve  4  pde's  using  Fourier  series 

C  Navier  Method! 

C 

C  iout  =  1  :  buckling  critical  stresses 

C  2  :  fixed  Ef,  h,  b/a,  varying  tf/tc,  Ef/Ec 

C  3  :  maximum  deflection  to  screen  only 

C  else  :  complete  deformation  output  at  (x,y,z)  &  p(x,y) 

C  


dimension  vcl(30,30,5),vc2(30,30,5),vc3(30,30,5),hfd(3) 

dimension  u3(30,30),s3(30,30),sl(30,30),s2(30,30),p(30,30) 

dimension  cm(4,4),bm(4),um(4),scfac(4),aa(4,5), 

+      pmn(19,19),u3mn(19,19),s3mn(19,19),slmn(19,19),s2mn(19,19) 

common  pi,a,b,hc,hf,vc,vf,ec,ef,po,rnx,nft 

open(unit  =  l,file  =  'seql.dat',status  =  'old') 

open(unit  =  2,file  =  'seql.out',status  =  'old') 

data  hfd/.25,.l,.05/ 

read(l,*)  a,b,hc,hf 

read(l,*)  vc,vf,ec,ef 

read(l,*)  iout,rnx,po,npx,npy,nft 

pi  =  3.14159 


C 


if(iout.ne.l)  goto  105 
write{2,101)  a,b 

101  format(5x,'Table  Generated  for  a=  ',f6.1,'  and  b=  ',f6.1/ 
&  /6x,'tf       tc',7x,'tf/b         tc/b',6x,'m     Scrl.min', 

&  2.\,'m     Scr2.min     Scr.min'/) 

do  104  i=l,3 
hf=hfd(i) 
do  102  j  =  1,11 
hc  =  real(j-l) 
call  screva 

102  continue 
writc(2,103) 

103  format(/) 

104  continue 
goto  999 

C — 

105  if(iout.ne.2)  goto  111 

hc=(10-4*hf)/2 
jcount=l 

106  ec  =  ef/(10*jcount) 
if(jcount.eq.3)  ec  =  ef/50 

111  do  120m  =  l,nft 
do  120n=l,nft 
pmn(ni,n)  =  0. 
u3mn(m,n)  =  0. 
s3mn(m,n)  =  0. 
slmn(m,n)  =  0. 
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s2mn(m,n)  =  0. 
120  continue 
C 

call  coefs(pmn,u3mn,s3mn,slmn,s2mn) 
do  200  i  =  l,npx+l 
xa=(i-l)/real(npx) 
do200j  =  l,npy+l 
yb  =  G-l)Aeal(npy) 
u3(i,j)  =  0. 
s3(i,j)=0. 
sl(i,j)  =  0. 
s2(i,j)  =  0. 
p(i,j)  =  0. 

do  140  m  =  l,nft,2 
col  =  sin(m*pi*xa) 
co2  =  cos(ni*pi*xa) 
do  140  n  =  l,nft,2 
co3  =  sin(n*pi*yb) 
co4  =  cos(n*pi*yb) 

u3(i,j)  =  u3(i,j)  +  u3mn(m,n)*col*co3 
s3(i,j)  =  s3(i,j)  +  s3mn(m,n)*col*co3 
sl(i,j)  =  sl(i,j)  +  slmn(m,n)*co2*co3 
s2(i,j)  =  s2(i,j)  +  s2mn(m,n)*col*co4 
P(i.j)  =  P(ij)  +  pmn(m,n)*col*co3 
140  continue 

do  150  k  =  1,3 
zh  =  (k-2)*hc 
vcl(i,j,k)  =  zh*sl(i,j) 
vc2(i,j,k)  =  zh*s2(i,j) 
vc3(i,j,k)  =  u3(i,j)  +  zh*s3(ij) 
150          continue 
200  continue 
C - 

if(iout.ne.2)  goto  230 

\vrite(2,210)  b/a.htyiiccf/cc 
write(2,220)  (100*vc3(i,2,2),i=  11,21) 
210      formal(5x,'b/a  =',f5.2,6x,'hf/hc  =',f6.2,6x,'Ef/Ec  =',f6.2) 
220      formal(lx,ll(f7.4)/) 
if(jcount.ge.3)  goto  999 
jcount=jcount  +  l 
goto  106 

C 

230  if(iout.eq.3)  goto  666 
do  300  i=l,npx+l 
do300j  =  l,npy+l 
write(2,600)  i,j,p(i,j) 

write(2,610)  (k,vcl(ij,k),vc2(i,j,k),vc3(i,j,k),k=  1,3) 
300  continue 
g^)to  999 
600  I"ormat(5x,2i6,5x,el0.2) 
610  format(3(i4,3el2.4/)) 
666  write(6,*)  (vc3(2,2,k),k  =  l,3) 
999  end 
C --- 


71 


subroutine  coefs(pmn,u3mn,s3mn,slmn,s2mn) 
dimension  cm(4,4),bm(4),um(4), 

+      pmn(19,19),u3mn(19,19),s3mn(19,19),slmn(19,19),s2mn(19,19) 
common  pi,a,b,hc,hf,vc,vf,ec,ef,po,rnx,nft 
gc  =  ec/2/(l  +  vc) 
gf=ef/2/(l  +  vf) 
dl  =  4*gf*hf**3/3/(l-vf) 
ac  =  2*(l-vc)/(l-2*vc) 
tl  =  8*gf*hf**2*hc/(l-vf) 
t2  =  2*gc*hc**2/3 
t3  =  hc*(l  +  gc*hc/(6*gf  *hf)) 
t4  =  hc*((l  +  vf)/(l-vf)  +  gc*hc*(ac-l)/(6*gf*hf)) 
t5  =  gc/(2*gf*hf) 
t6  =  2*hf/(l-vf) 
do  160m  =  l,nft,2 
ftl  =  m*pi/a 
do  160  n=l,nft,2 
do  110  i  =  1,4 
bm(i)  =  0. 
um(i)  =  0. 
do  110  j  =  1,4 
cm(i,j)  =  0. 
110  continue 

ft2=n*pi/b 
t7  =  ftl**2  +  ft2**2 

t8  =  ftl**4  +  2*(ftl**2)*(ft2**2)  +  ft2**4 
cm(l,l)=8*dl*t8  +  2*hc*gc*t7-2*rnx*ftl**2 
cm(l,3)  =  (2*hc*gc-tl*t7)*ftl 
cm(l,4)  =  (2*hc*gc-tl*t7)*ft2 

cm(2,2)  =  8*dl*hc*t8  +  t2*t7  +  2*ac*gc-2*rnx*hc*ftl**2 
cm(3,l)  =  ftl*(t6*t7-t5) 
cm(3,3)  =  -t3*t7-t4*ftl**2-t5 
cm(3,4)  =  -t4*ftl*ft2 
cm(4,l)  =  ft2*(t6*t7-t5) 
cm(4,3)  =  -t4*ftl*ft2 
cm(4,4)  =  -t3*t7-t4*ft2**2-t5 
pmn(m,n)  =  16*po/pi**2/m/n 
bm(l)  =  pmn(m,n) 
bm(2)  =  -pmn(m,n) 
if(hc.ne.O.)  goto  140 

u3mn(m,n)  =  bm(l)/cm(l,l) 
s3mn(m,n)  =  0. 
slmn(m,n)  =  0. 
s2mn(m,n)  =  0. 
goto  160 
140  call  simeq(4,cm,bm,um) 

u3mn(m,n)  =  um(l) 
s3mn(m,n)  =  um(2) 
slmn(m,n)  =  um(3) 
s2mn(m,n)  =  um(4) 
160  continue 
return 
end 

C      
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C  Subroutine  to  Evaluate  the  Critical  buckling  loads 

C 

subroutine  screva 

common  pi,a,b,hc,hf,vc,vf,ec,ef,po,riLX,nft 

gc  =  ec/2/(l  +  vc) 

gf=ef/2/(l  +  v0 

dl  =  4*gf*hf**3/3/(l-vf) 

ac=2*(l-vc)/(l-2*vc) 

tl  =  6*dl*hc/hf 

t2  =  2*gc*hc**2/3 

t3  =  hc*(l  +  gc*hc/(6*gf*hf)) 

t4  =  hc*((l  +  vf)/(l-vO  +  gc*hc*(ac-l)/(6*gf*hO) 
t5  =  gc/(2*gf*hf) 

t6  =  2*hf/(l-v0 
do  150  i=l,2 

m  =  l 

scrm  =  10.e+9 

101  ftl  =  m*pi/a 
ft2  =  pi/b 

t7  =  ftl**2  +  ft2**2 

t8  =  t7**2 

if(i.eq.2)  goto  102 

cll  =  8*dl*t8  +  2*hc*gc*t7 

cl3  =  (2*hc*gc-tl*t7)*ftl 

Cl4=(2*hc*gc-tl*t7)*ft2 

c31  =  ftl*(t6*t7-t5) 

c33  =  -t3*t7-t4*ftl**2-t5 

c34  =  -t4*ftl*ft2 

C41  =  ft2*(t6*t7-t5) 

c43  =  -t4*ftl*ft2 

c44  =  -t3*t7-t4*ft2**2-t5 

tmp  =  cl3*(c31*c44-c41*c34)-cl4*(c43*c31-c41*c33) 

scr  =  (cll-tmp/(c33*c44-c43*c34))/4/hf/ftl**2 

goto  103 

102  c22  =  8*dl*hc*t8  +  t2*t7  +  2*ac*gc 
if(hc.eq.O.)  goto  130 

scr  =  c22/4/hf/hc/ftl**2 

103  if(scr.ge.scrm)  goto  130 
scrm  =  scr 

m  =  m  +  l 

if(m.lt.50)  goto  101 
write(6,120)  i,m 
120      format(lx,'  In  SCR(',il,')  exceeded  max  m  =  ',i4) 

stop 
130      if(i.eq.l)  then 
scrml  =  scrm 
ml  =  m-l 
else 

scrm2  =  scrm 
m2  =  m-l 
endif 
150  continue 

scrmin  =  min(scrml,scrm2) 

write(2,200)  2*hf,2*hc,2*hf/b,2*hc/b,ml,scrml,m2,scrm2,scrmin 
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200  format(2x,2f7.2,lx,2fl0.5,2(2x,i2,2x,f9.0),2x,f8.0) 
return 
end 
C 

C  Subroutine  to  Solve  Simulteneous  Equations 

C 

subroutine  siineq(ne,cm,bm,um) 
dimension  cm(4,4),bm(4),um(4),scfac(4),aa(4,5) 
do  20  i=l,ne 
do  10  j  =  l,ne 
aa(i,j)=cm(i,j) 
10        continue 

aa(i,ne+l)  =  bm(i) 
20    continue 
do  50  i=l,ne 
big  =  abs(aa(i,l)) 
do  30  j  =  2,ne 

anext  =  abs(aa(i,j)) 
if(anext.gt.big)  big=anext 
30       continue 

if(big.lt..00QQ01)  then 
write(6,40)  i 
40  format(lhl,//,'**  Elements  in  Row  ',i2,'  are  Zeros  **') 

stop 
end  if 

scfac(i)  =  l./big 
50    continue 
do  60  i=l,ne 
do  60  j  =  l,ne+l 

aa(i,j)  =  aa(i,j)*scfac(i) 
60    continue 

do  110  i=l,ne-l 
ipvt  =  i 
ipl  =  i+l 
do  70  j  =  ipl,ne 

if(abs(aa(ipvt,i)).lt.abs(aa(i,j)))  ipvt  =  j 
70        continue 

if(abs(aa(ipvt,i)).lt..000001)  then 
print  *, 'solution  not  feasible.', 
&  '  a  near  zero  pivot  was  encountered' 

stop 
end  if 

if(ipvt.ne.i)  then 
do  80  jcol=l,ne+l 
temp  =  aa(i,jcol) 
aa(i,jcol)  =  aa(ipvt,jcol) 
aa(ipvt,jcol)  =  temp 
80  continue 

end  if 
do  100  jrow  =  ipl,ne 

if(abs(aa(jrow,i)).eq.O.)  goto  100 
ratio  =  aa(jrow,i)/aa(i,i) 
aa(irow,i)  =  ratio 
do  90  kcol  =  ipl,ne  +  l 
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aa(jrow,kcol)  =  aa(jrow,kcol)-ralio*aa(i,kcol) 
90  continue 

100      continue 
110  continue 

aa(ne,ne  + 1)  =  aa(ne,ne  +  l)/aa(nc,ne) 
do  130  j  =  2,ne 
nvbl  =  ne  +  l-j 
l  =  nvbl+l 

temp  =  aa(nvbl,ne  + 1) 
do  120  k  =  l,nc 

temp  =  temp-aa(nvbl,k)*aa(k,ne+l) 
120      continue 

aa(nvbl,ne  + 1)  =  temp/aa(nvbl,nvbl) 
130  continue 
do  140  i  =  l,ne 

um(i)  =  aa(i,ne+  1) 
140  continue 
return 
end 
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Listing  of  a  computer  program  written  in  FORTRAN  language  to  compute  and  print 
the  complete  defomation  data  of  a  simply  supported  rectangular  sandwich  plate 
under  uniformly  distributed  transverse  load,  based  on  Eringen's  theory. 

C  BENDING  OF  A  SIMPLY-SUPPORTED  RECTANGULAR  SANDWICH  PLATE 

C  UNDER  A  UMFORMLY-DIST.  TRANSVERSE  LOAD  ON  THE  UPPER  FACE 

C  AND  A  UNIFORMLY-DIST.  AXL\L  COMPRESSIVE  LOAD  ALONG  FACE 

C  EDGESX  =  0  ANDX  =  A. 

C 

DIMENSION  UP(30,30,5),VP(30,30,5),WP(30,30,5) 

DIMENSION  U(30,30),V(30,30),W1(30,30),W2(30,30),P(30,30) 

COMMON  PMN(19,19),UMN(19,19),VMN(19,19) 

COMMON  W1MN(19,19),W2MN(19,19) 

COMMON  PI,A,B,H,TC,TF,VC,VF,GC,GF,RIl,RIC,RIF,NFT,PO,RNX 
C       OPEN(  UNIT  =  1,  FILE  =  'ERINGEN1.DAT',  STATUS  =  'OLD'  ) 

OPEN(  UNIT  =  2,  FILE  =  'ERINGEN1.0UT',  STATUS  = 'OLD'  ) 

A  =100. 

B  =  100. 

TC  =  10. 

TF  =  2. 

VC  =  .3 

VF=.3 

EC  =  2000. 

EF=  10000000. 

RNX=10. 

PO  =  100. 

NFT  =  9 

NPX  =  4 

NPY  =  4 

PI  =  3.14159 
cc      jcount  =  l 
cc      READ(5,*)  TF 
cc      TC=10-2*TF 

C       READ(1,*)  A,B,TC,TF,VC,VF,EC,EF 
C       READ(1,*)  PO,RNX,NPX,NPY,NFT 
C 

cclO   ec  =  ef/(10*jcount) 
cc      if(jcount.eq.3)  ec  =  ef/50 

H=TC+TF 

RI1  =  TF**3/12. 

RIC  =  TC**3/12. 

RIF=TF*H**2/2. 

GC  =  EC/2./(l  +  VC) 

GF  =  EF/2./(l  +  VF) 

DO  100M  =  1,NFT 
DO  100  N=1,NFT 
PMN(M,N)=0. 
UMN(M,N)  =  0. 
VMN(M,N)=0. 
W1MN(M,N)=0. 
W2MN(M,N)  =  0. 
100   CONTINUE 

CALL  COEFS 
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c 

DO  200  I  =  1,NPX+1 
XA=(I-1)/REAL(NPX) 
DO200J  =  l,NPY  +  l 
YB  =  {J-1)/REAL(NPY) 
U(I,J)=0. 
V(I,J)  =  0. 
W1(I,J)=0. 
W2(I,J)=0. 
P(I,J)  =  0. 

DO  120M=l,NFr,2 
C01  =  SIN(M*PI*XA) 
C02  =  C0S(M*PI*XA) 
DO  120N  =  1,NFT,2 
C03  =  SIN(N*PI*YB) 
C04  =  C0S(N*PI*YB) 
U(I,J)  =  U(I,J)  +  UMN(M,N)*C02*C03 
V(I,J)  =  V(I,J)  +  VMN(M,N)*C01*C04 
Wl(I,J)  =  Wl(I,J)+WlMN(M,N)*C0rC03 
W2(I,J)  =  W2(I,J)+W2MN(M,N)*C01*C03 
P(I,J)  =  P(I,J)  +  PMN(M,N)*C01*C03 
120  CONTINUE 

DO  130  K=  1,3 
ZH  =  (K-2)*TC/2 
UP(I,J,K)  =  ZH*U(i,J) 
VP(I,J,K)  =  ZH*V(I,J) 
WP(I,J,K)  =  W1(I,I)  +  ZH*(2/h)*W2(I,J) 

130  CONTINUE 

200  CONTINUE 

cc      write(2,210)  b/a,tf/tc,cf/ec 
cc      write(2,220)  (100*wp(i,2,2),i=  11,21) 

CC210  format(5x,'b/a  =  ',f4.2,5x,'  tf/tc  =  ',f5.2,5x,'  Ef/Ec  =  ',15.2) 
cc220  format(lx,ll(f7.4)/) 
cc      if(jcount.ge.3)  goto  999 
cc      jcount=jcount  +  l 
cc      goto  10 
c       \vrite(6,*)  wp(2,2,2) 
DO  300  I  =  1,NPX  +  1 
DO  300  J  =  1,NPY  +  1 

WRITE(2,600)  I,J,P(I,J) 

WRITE(2,610)  (K,UP(I,J,K),VP(I,J,K),WP(I,J,K),K=  1,3) 

300  CONTINUE 
600  FORMAT(5X,2I6,5X,E10.2) 
610  FORMAT(3(I4,3E12.4/)) 
999  END 
C 

SUBROUTINE  COEFS 

COMMON  PMN(19,19),UMN(19,19),VMN(19,19),W1MN(19,19),W2MN(19,19) 

COMMON  PI,A,B,H,TC,TF,VC,VF,GC,GF,RIl,RlC,RIF,NFT,PO,RNX 

TO  =  RNX/GC 

T1  =  GC*RIC/GF/RIF 

T11  =  RIC*TF/RIF/TC 

T2  =  2*(1-VC)/(1-2*VC) 
T3  =  (1-VF) 
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T33  =  GF*RI1/T3/GC/RIC 
DO  100  M=1,NFT,2 

T4  =  (M*PI/A) 

T5  =  (T4**2)*T0*TC 

DO  100  N=l,NFr,2 

T6  =  (N*PI/B) 

T7  =  (T4**2  +  T6**2)*TC**2 

T8  =  (2/T3  +  T1*T2)  *T7/12.  +  Tl 

DIMN  =  T7*(-T1  +  (T33*T7/3.  +  l.)*T8)-2.*T5*T8 

D2MN  =  T7*(T33*T7/6.  +  Tll)  +  12.*T2*Tll-T5 

PMN(M,N)  =  16.*PO/PI**2/M/N 

UMN(M,N)  =  -(PMN(M,N)/2./GF/DlMN)*(RIC/RIF)*(M*PI*TC/A)*H 

VMN(M,N)  =  UMN(M,N)*A/B 

W1MN(M,N)  =  (PMN(M,N)*TC/GC/D1MN)*T8 

W2MN(M,N)  =  PMN(M,N)*TC/2./GC/D2MN 
100   CONTINUE 
RETURN 
END 
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Listing  of  a  computer  program  written  in  FORTRAN  language  to  evaluate  the 
critical  stresses  for  instability  of  a  simply  supported  rectangular  sandwich  plate,  based 


on  Eringen's  theory 


DIMENSION  TFD(3) 

COMMON  A,B,VC,VF,EC,EF 

DATA  TFD/.5,.2,.l/ 

VC  =  .3 

VF=.3 

EC  =  2000. 

EF=  10000000. 

A  =100. 

B  =  100. 

OPEN(UNIT  =  3,FILE  =  'ERIN-1.0UT',STATUS  =  'UNKNOWN') 

WRITE(3,101)  A,B 

101  FORMAT(5X,'TabIe  Generated  for  a=  ',F6.1,'  and  b=  ',f6.1/ 
&  /6X,'tf       tc',7X,'tf/b         tc/b',6X,'m     Scrl.min', 

&  2X,'m     Scr2.mm     Scr.min'/) 

DO  104  1  =  1,3 
TF  =  TFD(I) 
DO  102  J  =  1,11 
TC  =  2*REAL(J-1) 
CALL  SCREVA(TC,TF) 

102  CONTINUE 
WRITE(3,103) 

103  FORMAT(/) 

104  CONTINUE 
STOP 

END 


CC 


SUBROUTINE  SCREVA(TC,TF) 
COMMON  A,B,VC,VF,EC,EF 
PI  =  3.14159 
GC  =  EC/2/(l  +  VC) 

GF  =  EF/2/(l  +  VF) 

C1  =  (PI**2*EF)/(12*(1-VF**2)) 

C2  =  (1-VF)/(PI**2)*(GC/GF) 

C3  =  (1-VF)*(2*(1-VC)/(1-2*VC))/12*(GC/GF) 

C4=  12*(2*(1-VC)/(1-2*VC))/PI**2 

DO  150  1  =  1,2 

M  =  l 

SCRM  =  10.e  +  9 

101  X1  =  A/B/M 
X2=TC/B 
X3=TF/B 

SCRF  =  C1*(X3**2)*(1/X1  +  X1)**2 
IF(I.EQ.2)  GOTO  102 

ETA1  =  (3*C2*X2/X3**3)/(1  +  1/X1**2+C2*X2/(X3*(X3  +  X2)* 
&  +C3*X2**3)) 

SCR  =  SCRF*(1  +  ETA1) 
GOTO  103 

102  IF(TC.EQ.O.)  GOTO  130 

ETA2  =  C2*(X2/X3)**3*(1  +  X1**2  +  C4*(X1/X2)**2)/ 


79 


&  (1/X1  +  X1)**2/(X2  +  X3)**2 

SCR  =  SCRF*(1  +  ETA2) 
103      IF(SCR.GE.SCRM)  GOTO  130 
SCRM  =  SCR 
M  =  M  +  1 

IF(M.U.50)  GOTO  101 
WRITE(6,120)  I,M 
120      FORMAT(lX,'  In  SCR(',I1,')  exceeded  max  m  =  ',14) 

STOP 
130      IF(I.EQ.l)  THEN 
SCRM1  =  SCRM 
M1  =  M-1 
ELSE 

SCRM2  =  SCRM 
M2  =  M-1 
ENDIF 
150  CONTINUE 

SCRMIN  =  MIN(SCRM1,SCRM2) 

WRITE(3,200)  TF,TC,X3,X2,M1,SCRM1,M2,SCRM2,SCRMIN 
200FORMAT(2X,2F7.2,1X,2F10.5,2(2X,I2,2X,F9.0),2X,F8.0) 
RETURN 
END 
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